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SECTION 1.0 
INTRODUCTION 


The primary objective in the development of the SEAHT 
program is the improvement of sea surface height values obtained 
from altimeter data along intersecting arcs (i.e., arcs whose 
groundt racks cross) . The basic concept in the improvement 
scheme is that, at the crossovers, the altimeter measurements 
are to the same point on the surface and crossover differences 
are due to differences in: 

• orbit error 

• tides 

• currents 

• propagation effects 

• altimeter noise and bias 

I 

• eddies and other temporal sea surface height 
effects 

Of these effects, all except orbit error can be corrected or 
have amplitudes below the approximately 1 ra level, and the 
effects of measurement noise and propagation effects are expected 
to be reduced to near the 10 era level. The satellite orbital 
height error is expected to be the dominant error source in 
most cases. Over short time periods (less than a quarter of 
a revolution) , the orbit error can be modeled as a polynomial 
in time, and a set of orbit error parameters can be estimated 
such that the set of intersections are as consistant as allowed 
by unmodeled (or improperly modeled) temporal sea surface height 



variations. In general, this orbit adjustment would be ex- 
pected to be made relative to one or more passes whose orbits 
were known to be very good. This allows the estimated sea 
surface heights to be given an ,ahsqlut,e- interpretation rather 
than being only a set of relative heights. 

Given a set of arcs, algorithms have been deyeloped to 
categorize and find the intersections of altiijieter passes. 
Algorithms are also developed for the estimation of the orbit 
error parameters using a Bayesian least squares technique. 


2 



SECTION 2 . 0 

CROSSING ARC -ADJUSTMENT ANALYSIS 


2.1 FORMULATION OF CROSSOVER DIFFERENCES 

In the study of sets of intersecting arcs, it is necessary 

to construct matrices relating to a pass at the time that it 

crosses the groundtrack for another pass. This will be accomplished 

through double subscripts, with the first denoting the pass of 

interest, and the second the pass which is being intersected. 

E.g., gi i then refers to the i^^ pass at the time of intersection 
^ J 

of the groundtrack from the j pass. In all cases that will 
be considered, the arc lengths will be chosen sufficiently short 
so that the intersections are unique,. and one arc does not cross 
another more than once. 

( 

Assuming that tide corrections have been made to altimeter 
sea surface height measurements, the result is the altimeter 
geoid measurement contaminated by, in general, orbit error. 

This geoid height measurement can be expressed as 



g 


ij 


Ag. . 


+ n.. 


where 


CD 


is the geoid height measurement 

g^. is the actual geoid height 

1 

Agij is the error in the geoid height due to satellite 
height error plus any biases in the altimeter 
measurement 
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it may 


where 


point ( 


3ut , by 
ffe have 


rhere d 
It the 


n. . is the error in the measurement due to noise 
in the altimeter measurement and to temporal 
sea surface height effects. 


Although no analytical model is available for, g^ 
be represented empirically by a polynomial in time. 


' ®i0 '"ij * ^12 + 


- K. 


C2) 


t^j is the time along the i^^ arc at point (i,j) 

Considering now the geoid height measurement at the 

th 

i,j) along the j arc, the model equation is 
g-- = g* • + As-~ + n-. 

definition g^j* = Therefore, using (1) and (3), 


- 


®ij 


“ g • • 


= Ag . . - 

6, 1 


Ag . . 


+ n. 


11 


- n 


11 


j^j is the difference between the geoid height measurements 
Ci,j) point. Substituting from Eqn. (2), 


^i 


K. 

1 


. . = > a. tf . - y a.^ tf, + n. . - n.. 
11 A-' 1® 11 ^ 1® 11 11 11 

e=0 e=0 


(4) 
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Since Equation (4) is expressed in terms of geoid height 
differences, a constraint must be used to anchor the solution. 
For a one parameter solution, one constrained revolution is 
sufficient. All other revolutions are adjusted such that the 
differences in heights between each of these and the constrained 
revolution will be minimized. Since orbital errors have 'pre- 
dominantly a once per revolution type of variation, a low 
order polynomial should serve as a satisfactory model for short 
arcs. In particular, a quadratic would be expected to be the 
maximum required order for arc lengths ’not exceeding’ a quarter 
of a revolution C~25 minutes for GEOS- 3) and shorter arcs should 
require only a first (or zero) 'order polynomial. 


The objective is then to solve for the set of coefficients 
which will minimize the in some sense, subject to the 

constraint that certain arcs are not to be adjusted. The con- 
straints can be applied through several different techniques. 

We choose the- Bayesian approach since it is consistent with 
available software and amounts to nothing more than supplying 
additional measurements to the estimation scheme. The least 
squares technique will' be chosen for the coefficient estimation. 
Th4.s procedure will also- be minimum variance if the noise' 
difference, n^^ - is uncorrelated from one crossing to 

another and its variance is used to weight the intersection 
data. The degree to which the errors are uncorrelated depends 
primarily upon the correlation between temporal sea surface 
height variations but, in any case, correlation information 
is not expected to be available. The least squares technique 
then requires' the minimization of 


N N 


s = 2) 2 

i“l j“i+l 



K-1 


le 


e=0 




C5) 




s 



where 


N is the total iiximber of arcs 

= 1 if arc j intersects i 
= 0 if not 


w 


• .= weighting factor based on n. . - n.. 

to t. 

- l/E [n.. - n..]2 


with E the expected value operator. Additionally, we have 
assumed that 


so that the number of orbit error parameters is taken to be 
the same for all arcs. This is not an essential assumption, 
but substantially simplifies the software implementation. 
Constraints can be used to effectively vary the number of 
adjusted parameters from one arc to the next. 


2.2 IMPLEMENTATION CONSIDERATIONS 

Before discussing the solution of C5) it is, pertinent 
to consider the number of intersections required to generate 
a solution of a given order. First, note that the total 
number of intersections is 

f N N 

"i = E E hr 

i=l j=i+l 
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whereas the number of coefficients to be determined is 
N 

^ K = NK (7) 

. i=l 

t 

Obviously, for this method to be practical we must require 
that 


Ni > 


or 


N N 


L E 

i=l j=i+l 


6 . . 

\3 


> NK 


( 8 ) 


In addition, we cannot solve for K coefficients along 
the i^^ arc if there are fewer than K arcs intersecting arc 
i. Therefore, we also require 


N 

i=i 

In order to assess the number of arcs which it may be 
feasible to process 'in a computer, consider the case of two 
sets of arcs with each arc in the first set intersecting each 
arc in the second set. We will denote the number of arcs in 
the first set by N and the number of' arcs in the second set 
by N^, For a solution, we require that the number of inter- 
sections be greater than or equal to the number of coefficients , 
or 
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from which 



-N^-K 


( 10 ) 


Minimum values of N for various values of K and Ni are listed 

a D 

in Table 1. 

Minimum 

^ ^a 

12 2 

13 2 

15 2 


2 

2 

2 

2 

2 


3 

4 

5 

10 

20 


6 

4 

4 

3 

3 


3 4 12 

3 5 8 

3 10 5 

3 20 4 

3 50 4 


Table 1. Minimum Numbers of Crossing Arcs 
for Different Numbers of Adjusted Parameters 

Note that the limiting value of N for large values of N, is 

D 

K+l. This can be deduced from Equation (10), since 


£im 

Nb^» 


Nb-K 


K 


( 11 ) 


For all finite values of N|^, the ratio is greater than K so 
No must be at least K+l, 

Ci 
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Looking now at the approximate core requirements, we 
consider the normal matrix to be the array which primarily 
limits the number of arcs' which can be simultaneously processed. 
Allowing for this array to be double precision, but with only 
the diagonal and above elements stored, the niomber of words 
required for the M dimensional normal matrix is 

MCM+1) 

N^, = (2) = M(M+1) 

^ 2 

Substituting M = K(N^ + N|^) for the number of parameters, the 
core required is 

Nt = K^CN^ + +- K(N^ + -(12) 

Table 2 shows the core required for the normal matrix for two 
different relations of to N. . In one case, we set N = N, , 
and in the other case we set = 1 as is appropriate for the 
program implementation in which the arcs in one direction are 
partitioned, as will be discussed in Section 2,3, 

Considering the maximum core available on the Wallops 
computers to be on the order of 70K words, it is evident from 
Table 2 that the limitation on for the unpartitioned solution 
is approximately 100 passes for a one parameter solution and 
approximately 50 passes' for a two parameter solution. For a 
partitioned solution, these limits are almost doubled, allowing 
approximately 200 one parameter arcs or 100 two parameter arcs. 
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K 

N . . 

a 

.. N 

T 


t • » 

■ ' ,Nb 

H 

II 

1 

2 

20 

12 


5 

110 

42 


10. 

420 

132 


20 

1640 

462 


50 

10100 ■ 

2652 


100 

40200 

10302 


150 

90300 

22952 


200 

160400 

40602 

2 

4 

2 72 

110 


10.. 

1640 

506 


20 

6480 

1806 


50 

40200 

10506 


100 , 

- 1 

160400 

41006 

3 

6 

1332 

462 


10 

3660 

1122 

. - 

20 

14520 

4032 


50. , 

90300 

23562 


100 

360600 

^ , 92112 


f . t 


Table 2. Normal Matrix Storage Requirements 
for Different Numbers of Arcs and Parameters Per Arc 
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2.3 


LEAST SQUARES ESTIMATION OF ORBIT IMPROVEMENT 
PARAMETERS 


The least squares estimate of the parameters which 
minimize the altimeter crossovers is obtained by finding the 
set of a’s which minimize the sum of squares of crossover 
differences given by Equation (5). Equivalently, we can con- 
sider the least squares solution of the set of linear equations 
given by Equation (4) , which may be written symbolically in 
matrix form as 


7 = BT + e -C13) 

where 

Y is the vector of geoid height differences. 

B is the matrix of partial derivatives of the geoid 
height differences with respect to the orbit error 
parameters. This matrix will consist of powers of 
time from some reference time for each pass. 

Y is the vector of orbit error coefficients to be 
solved for. 

Y is the vector of errors in the crossover differences 
due to measurement noise. 


Table 3 shows Equation (13) written out explicitly for a two 
coefficient parameterization of orbit error and 4 North-South 
passes crossing 4 South-North passes.* For convenience, the 


*The solution of the set of equations given in Table 3 is 
actually singular. Intuitively, it is obvious that whatever 
solution is obtained for the ajo’s, the same constant added 
to each of these coefficients is an equally valid solution. 

The solution can be made non-singular through the use of a 
constraint on one or more of the passes. We will assume that 
such information is to be used, but will not explicitly include 
it in the equations shown here, 

11 omginal 
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Table 3 . Matrix Form o£ Linear Equations to be Solved for 
4 Passes Crossing 4 Passes in the Opposite Direction and 

2 Orbit Parameters Per Pass 
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North- South arcs have been denoted as arcs 1-4 and the South- 
North passes as arcs 5-8. 

The weighted least squares solution of Equation (13) 

is [1] 


Z = b'^WY 


( 14 ) 


where 


A 

Z is the estimated set of coefficients 

W ^ = E[e e^] = covariance of the measurement noise, 
with E the expectation operator. 

That the solution allows partitioning can be seen by 
rewriting the equations given in Table 3 in the form: 



13 
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Assuming only that the n^'s are uncorrelated, and denoting 


w:^ = E[n^ nT], 


C16) 


Equation (14) for the 8 arc 16 parameter case of Table 3 becomes 


^1 



0 

0 

0 

B>iBl0 

/S 

^2 
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0 
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®2'*2®20 

/N 

^3 

= 

0 ■ 

0 

®3^3^3 

0 

B3**3®30 

A 

^4 


0 

0 

0 

®>4®4 

®I'*4®40 

A 

_^0_. 




(B^W3B3o)T 


1=1 


r 

1 





X 


b"^w y 
2 2^2 

T 

B W Y 
^3'^3 3 

bIV4 




(17) 


} 

The matrix inverse in Equation (17) can be obtained via parti- 
tioning [1] with the resulting expressions for the coefficient 
sets found to be 


^0 ^0 


^[bToWiY. - bT^w.b.cbTw.b.)-' bTw.y,] 

i=l 


( 18 ) 
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(19) 


ai = (Bjw.Bp-l B^W-Y. - (B^W-Ep-l Bjw.B.„ao 


where 


M. 


*T 

= E 




i=l 


The extension of- the solution given by Equations (18-20) 
to any number of arcs is achieved through two directions. Let 
us refer to the arcs for which the a^ parameters have been 
assigned as North-South arcs, and the set of arcs for which 
the Eq set of parameters have been assigned as South-North 
arcs. To increase the number o£ North- South arcs, we simply 
increase the summation limit in Equations (18) and (10), To 
increase the number of South-North arcs, we must increase the 
number of Eg parameters and thus the size of Mg which is, in 
general, the matrix providing the major limitation on the 
number of passes which can be processed in one computer run. 

The estimation scheme implemented in SEAHT is the par- 
titioned least squares type solution described above, with 
the number of adjustable parameters per pass selectable as 
1, 2, or 3, North-South arcs are partitioned out', so virtually 
any number of North-South passes can be processed simultaneously. 
South-North arcs provide the "common" parameter set, and the 
number of such passes is limited by available computer core. 
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SECTION 3.0 

PROGRAM IMPLEMENTATION 


3.1 PROGRAM ORGANIZATION 

The SEAHT program utilizes overlays, and consists 
primarily of a main program (plus common block initialization) , 
an input set of routines, and an estimation set of routines. 

The last set of routines, however, is divided into two groups, 
and there is an additional segment used for core allocation. 
Thus, there are a total of four program "LINKS" plus a main 
program segment. 

The main program is responsible for the flow of the 
program, and for the computation of necessary array sizes . 

In LINK 1, optional input is read from cards and input data 
is read from tape(s), and the arc crossings are computed. 

The crossing data is stored on disk in the form to be used by 
the estimation routines. The order of the crossings is: all 
crossings of the first North- South arc, followed by all cross- 
ings of the second North-South arc, etc. This format allows 
the partitioning scheme discussed above to be implemented. 

LINK 1 contains routines with fixed array sizes, with dimensions 
set to accommodate the maximum amounts of data per arc and 
maximum number of arcs that are expected. However, these 
dimensions do not result in excessive core utilization. 

At the conclusion of the computation of the arc cross- 
ings, the needed sizes for all arrays in the estimation process 
are calculated by the main program. LINK 1 is then overlaid 
by LINK 2 which contains a portion of the estimation routines. 
LINK 3 is then called to perform dynamic memory allocation 
for all arrays needed in the parameter estimation process. 

After this allocation, LINK 3 is overlaid by LINK 4 containing 


ORIGINAL 
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the remainder o£ the estimation routines. No further overlays 
are performed. 

Control of the estimation process is from subroutine 
SOLVE, which is responsible for controlling the normal matrix 
set-up, the application of constraints, the parameter adjust- 
ments, and the application of the estimated parameters to the 
input data, with the corrected, data then saved on tape. 

3.2 PROGRAM INPUT 

Input to the SEAHT program consists of one or more 
data tapes, and card input which requests certain program 
options. The latter category consists of a set of option 
cards, at least one of which is required. 

3.2.1 Data Tape Input 

Data input files to SEAHT consist of a header record 
followed by a set of data records, as specified on the follow- 
ing two pages. The header record identifies the time of the 
data by both the revolution (and segment) number and date. 

At the present time, the date is not used, while the revolution 
and segments numbers are used for the selection and deletion 
of particular data passes, as discussed below. 

Input tapes may be prepared in several different ways. 
They may be prepared by extracting information contained on 
each frame of an ARC [2] tape, in which case sea surface heights 
are major frame i'^2.0 or 3.2 seconds) averages. In addition, 
the ALTKAL program [3] produces a tape which is directly 
readable by SEAHT. This tape contains smoothed sea surface 
heights at the 1/second data rate. 
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In the current program implementation, only Words 1-5 
o£ the header record and Words 1-4 o£ the data record are read 
by SEAHT, and the £irst 3 words o£ the header record are not 
used. Words 4-5 o£ the header record and the £irst 4 words ' 
o£ each data record are critical to program operation. 
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SEAHT Input Format 
Header Record 


Word No. 

Type 

Description 

1 

I 

Year and Day of year (YYDDD) 

2-3 

DP 

Time of day (seconds) for first 
data record 

4 

I 

Segment number for altimeter data 
pass 

5 

I 

Rev number for altimeter data pass 

6 

I 

Measurement type 

40 = Global mode 

41 = Intensive mode 

7-19 

I 

Zeroes (used for making all records 
same length) 
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SEAHT Input Format 
Data Record 


Word No. 


1 

2 

3 

4 

5-13 


14 

15 


16 

17-19 


Type 

R 

R(LAT) 

R(LON) 

R 

R 


R(THITB) 
R(TREF) ^ 


R(RAGCAV) 

R 


Description 

Time from first data point (minutes) 

Latitude of sea surface height 

Longitude of sea surface height 

Sea surface height (meters) 

Not used. (For tapes created by 
the ALTKAL smoother, words 4-13 
contain sea surface heights at the 
10/sec data rate.) 

Tide height (meters) 

Tropospheric refraction correction 
(meters) 

Average AGC 

Zeroes for future additions 


Note: The *’sea surface height” valuers) have been corrected 

using the tropospheric refraction correction and tide 
height on the data record. Strictly speaking, the sea 
surface height is the altimeter measurement of geoid 
height, although including the effects of currents and 
other non- stationary sea surface height effects. 
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3.2.2 Option Card Input 


Through the use o£ one or more o£ nine di££erent keyword 
option cards, various aspects o£ the operation of SEAHT can be 
altered from the default mode. These option cards and the 
purpose for which they are nitilized are 

• AREA - defines the geographic area for which arc 

crossings are to be selected. 

• PARAM - selects the number of coefficients to be 

used for each pass to parametrize orbit 
error. 

• NODATA - specifies that arc crossings have been 

saved from a previous program execution 
and are to be input in place of the 
normal data tapes. 

• CONSTR - specifies that a particular pass is to be 

treated as having a good orbit with the 
parameters constrained in the adjustment 
process. 

• EDIT - specifies the level at which arc crossing 

differences are to be edited from the 
parameter estimation. 

• SELECT - specifies that only particular passes 

are to be selected from a' particular 
data 'tape. 

• DELETE -' specifies that all passes are to be 

selecteci from a particular tape except 
for those which are designated. 
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NAME 

COLUMNS 

FORMAT 

DESCRIPTION 


AREA 

1-4 

/ 

A4 

Defines the area of interest 


11-25 

E15.5 

Minimum acceptable value 
latitude 

of 


26-40 

E15.5 

Maximum acceptable value 
latitude 

of 


41-55 

E15.5 

Minimum acceptable value 
longitude 

of 


56-70 

E15.5 

Maximum acceptable value 

of 


longitude 

IF CARD OMMITTED; Default values are: minimum latitude, -90.0 degrees; 

maximum latitude, +90.0 degrees; minimum 
longitude, 0.0 degrees; maximum longitude, 
original page is -360.0 degrees. 

OR POOR QUALITY 
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PARAM 
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MfCC-SOSB 


NAME 

COLUMNS 

FORMAT 

DESCRIPTION 

PARAM 

1-5 

AS 

Indicates the number of 
parameters which are to be 
used for each pass. 


7-10 

14 

The number of adjusted 
parameters for each pass. 
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WFrr-5MI 


NAME 

COLUMNS 

FORMAT 

DESCRIPTION 


NODATA 

1-6 

A6 

Indicates crossing data 
exists 

already 


7-10 

14 

The logical unit number 

o£ the 


disk or tape where the crossing 
data is stored 


ORIGINAL RAGE IS 
OF POOR QUALITY 
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CONSTR 


CDNSTRI 
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■ I ] I i I 
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1 

y 


NAME COLUMNS FORMAT 

CONSTR 1-6 A6 


7-10 14 


11-25 E15.6 


DESCRIPTION 

Constrains selected revs (other 
revs are adjusted to fit 
constrained rev and constrained 
revs are not adjusted) 

Type of constraint 
= 1, pass totally constrained 

(bias, slope and curvature) 

= 2, bias adjusted, slope and curva 
ture constrained 

= 3, bias and slope adjusted, curva 
ture constrained 

Number of a revolution to be 
constrained 
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CONSTR 

(Cont.) 


NAME 

COLUMNS 

FORMAT 

DESCRIPTION 



CONSTR 
(Cont , ) 

26-40 

E15.6 

Segment number 
the rev number 

associated 
in columns 

with 

11-55 


Note: The use of a constraint is meaningless unless the parameter 

has been otherwise requested to be adjusted. B.g., the 
constraint of curvature for a pass is meaningless if only 
bias and slope are being adjusted, and should therefore 
not be requested. 
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EDIT 



NAME 

COLUMNS 

FORMAT 

DESCRIPTION 

EDIT 

1-4 

A4 

Sets editing criteria. 


11-25 

E15.6 

Maximum acceptable difference 
in sea surface heights at 
any crossing point (meters). 

IF CARD 

OMITTED: 

Default value 

is 100.0 meters. 
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DATA 


/ 


DATA 
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name columns format 


DESCRIPTION 


DATA 1‘4 A4 


Indicates the source of the data to 
be read in 


7-10 14 


Logical unit from which data 
is to be read 


11-25 E15.6 Number of files to be read 

from the logical unit designated 
in columns 7-10 


As many as six of these cards are possible barring limitations 
set outside the program. 
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SELECT 


SELCCTl 22} 4763. 
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NAME COLUMNS FORMAT 

SELECT 1-6 A6 

7-10 14 

11-25 El 5. 6 

I 

26-i40 E15.6 

41-55 E15.6 


DESCRIPTION 


Selects revolutions from 
a specific tape 

Logical unit number of the 
tape from which selection is 
to be made 

Revolution to be selected 

Segment number associated with 
revolution designated in 
’ Columns" 11-25 

Revolution to be selected* 


* - if "not more than one revolution is selected, columns 

*1 j • > 

'41'- 8 0 should be blank. 
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SELECT (cont.) 


NAME COLUMN FORMAT 


DESCRIPTION 


56-70 E15.6 Segment number associated 

with the revolution designated 
in columns 41-55. 


LIMITATION: Multiple SELECT cards are permitted, but no 

more than 50 selections may be made from' the 
same logical unit number. 


OF POOB 


31 


DELETE 



NAME COLUMNS FORMAT DESCRIPTION 

DELETE 1-6 A6 Deletes revolutions not desired 

from a specific tape 

7-10 14 Logical unit number of the 

tape from which deletion is 
to be made 

11-25 E15.6 Revolution to be deleted 

26-40 E15.6 Segment number associated with 

the revolution in columns 

41-55 E15.6 Revolution to be deleted* 

* If not more than one revolution is deleted columns 41-80 
should be blank. 
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DELETE 


NAME COLUMN FORMAT 


DESCRIPTION 


56-70 E15.6 Segment number associated 

with the revolution designated 
in columns 41-55 


LIMITATION: Multiple DELETE cards are permitted, but' 

no more than 50 deletions may be made from 
the same logical unit number. 
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ZERO 
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NAME 

COLUMNS 

FORMAT 

DESCRIPTION 

ZERO 

1-4 

A4 

Sets a longitude different than 
Greenwich as the longitude 
origin. 


11-25 

E15.6 

Longitude, relative to Greenwich 
to be used as origin. 

PURPOSE : 

Since the 

Tiiethod 

t w 

of finding intersections of* 


revolutions is dependent upon a decreasing longitude 
in the area of interest, it may be necessary to, 
temporarily, zero-set a certain longitude in the area 
of least interest. This is done within the program 
and requires no additional computation by the user. 
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MINUTE 
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MS.-C-5Q41 - 


NAME COLUMNS FORMAT DESCRIPTION 

MINUTE 1-6 A6 Permits the time to be input 

in minutes 


IF CARD OMITTED: The units of time are assumed to be in 

days and are changed to minutes within the 
program. 


ORIGINAL 
OF POO^ QUALPTY 


35 



3.3 


DATA SELECTION AND PARAMETER ESTIMATION PROCESS 


The estimation process begins with the acceptance o£ 
data that was collected in the area of interest. The revolu- 
tions are classified as being north-south or south-north. A ' 
revolution whose latitude increases to .the maximum possible 
value (approximately, 65.°) and, then, decreases is separated 
into a south-north segment and a north-south segment. This 
also applies to revolutions whose latitude decreases to the 
minimum value (approximately -65.°) and, then, increases.- 

Preliminary to -fin'ding crossing points, tests are made 
based on the fact that there will not be an intersection if 
the exiting latitude of a- south-north revolution is' less than 
the exiting latitude of' the north- south revolution or if the 
entering latitude of the north- south revolution is less than 
the entering latitude of the south-north revolution. If the 
exiting longitude of the -north'TSOUth revolution is_ greater 
than the entering longitude of the south-north revolution or 
the exiting longitude' of the south-north revolution is greater 
, than the entering longitude of the north-south revolution, 
there will not be a crossing. 

Every s outH- north' ‘revoflut ion is compared to every north 

* V ^ • - * 

south revolution to determine crossing point's by the followxng 
method: 

1. Equations for lines connecting the entering and 

exiting points of the revolutions are determined, 
and the intersections of these lines computed. 
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2 . 


The latitudes , corresponding to the longitude at 
the intersection, are found for each revolution, 
and the difference between these values is computed. 
If the absolute value of this difference is less 
than .0001 degrees, the average of the latitudes is 
taken, and a crossing point is declared at this 
average latitude and the longitude corresponding 
to the intersection point of the lines connecting 
entering and exiting points of the revolutions. 

/ 

3. If the difference is too great and the latitude 
is larger for the north-south revolution, these 
points become the entering points of the revolutions. 
Otherwise, they become the exiting points. This 
process continues until a crossing point is found 
or the number of iterations exceeds ten. 

The times into the revolutions and the sea surface heights 
at these crossing points are found by interpolation. 

A process to determine whether each revolution has crossed 
other revolutions a sufficient number of times is performed. 

The process repeatedly deletes revolutions and checks for 
sufficient crossings until it is not necessary to delete any 
additional revolutions . 


Examining one crossing point at a time, the difference 

in sea surface heights (sea surface height of the north-south 

revolution minus the sea surface height of the south-north 

revolution) at the crossing point is computed as well as the 

powers of time along the revolutions at this crossing point. 

T 

The elements of the normal 'matrix B WB are computed and placed 

T 

in their appropriate position as are the elements of the B WY 
matrix. The noise on each revolution is assumed to be the same. 


ORIGINAL 
OF POOR 
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Next, the no±.th- south revolutions* which -are to be con- 

T 

strained are noted and e'lemen'ts of the B WB matrix are altered 
to ensure that the 'coe£-£icien'ts of the bras- for this revolution 
will be small and that adg us tments on other revolutions will, 
necessarily, be made relative to.'this’ Tevolution. 

Handling all .crossing point, data for one north-south 
revolution at a time, inversion of the normal matrix is performed 
and biases are computed for the north-south revolutions based on 
the information from each particular north-south revolution under 
consideration with the south-north passes held fixed,. After 
all the north-south passes have been processed, the south- 
north constraints (if any) are applied and the normal matrix 
for the southTnorth passes is calculated (Eqn. 2 0)-. The south- 
north biases are then computed (Eqn. 18) and the north-south 
biases updated (Eqn. 19) . 

All data are corrected for the biases found. • 
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3.4 


PROGRAM OPERATION 


A source deck and binary decks of SEAHT are located 
at the Wolf Research and Development Group office in Pocomoke 
City, Maryland. The Honeywell 625/635 run setup consists of 
the appropriate control cards, the program on binary cards, 
the desired data cards, and a source of data input. The 
standard setup for a run consists of the following: 

• Specification of the FORTRAN FILE units for both 
the disk files used by the program and the tapes 
used for data input. 

• Program binary deck. 

® The set of option cards required for the particular 
run. 

A sample job setup is given in Section 3.4.2. 

3.4.1 FORTRAN Logical Units 

In addition to the normal card input and printer output 
files, SEAHT uses 5 scratch files, writes one output tape, 
writes 2 additional printer files, and normally uses 1-6 input 
tapes. All scratch files are sequentially accessed and logically 
could also be assigned to tapes. However, some of the files - 
particularly unit 12 - are rewound a number of times and could 
lead to inefficient program operation if assigned to tape. 

The FORTRAN logical units used are as follows; 


ORIGINAL PAGE IS 
OP POOR QUAury 
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LOGICAL USE 

UNIT 


5 

6 

10 

11 

12 

13 

14 

15 
18 

19 

* 


Option card input 
Standard printer output 
Output tape of corrected data 
Dis]c storage during execution 
Dislc storage during execution 
Dislc storage during execution 
Disk storage during execution 
Disk storage during execution 
Printed output 
Printed output 

Input tapes (normally multi-file) 


*Tliis set must be comprised of unique numbers different from 
the other units listed above. Multi- file tapes must be non- 
standard label. Single file tapes may be either standard 
label or non-standard label. A total of up to 6 input tapes 
may be used in a job. 

3.4.2 Sample Deck Setup 

The following is an example of a working deck setup. 
This particular setup would require that nineteen files be 
read from the tape with logical unit number 21 (in this case, 
tape #4572). A height bias would be found for all passes 
within the area: latitude; 45°-60®, longitude; 285°-300 , 

which are found to cross another pass. Rev 4476, segment 169 
would not be adjusted and the biases determined for the other 
passes would be relative to this constrained pass. 
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COL 


1 

$ 

$ 


$ 


$ 


$ 


$ 


8 

I DENT 
OPTION 


LINK 


LINK 


LINK 


LINK 


16 

112293, SEAHT 
FORTRAN 


LINKl 


LINK2,LINK1 

LINKS 


LINK4, LINKS 



PROGRAM 

BINARY 

DECK 


$ 

$ 

$ 

$ 

$ 

$ 


EXECUTE DUMP 

LIMITS 20,47K,10K 

TAPE 9 10 ,X6D, , , ,WJTAPE 

FF I LE 10, NSTDLB , NOSRLS , F I XLNG/ S , BUFS 1 2/ 5 , MLTFIL 

TAPE 21,X2D,,4572 

FF I LE 21, NSTDLB , NOSRLS , F I XLNG/ 1 0 , BUFS 1 Z/ 3 2 0 , MLTF I L 
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COL. 


1 

8 

16 

30 

45 


$ 

FILE 

11, DIR, SOS 




$ 

FILE 

12,D2R,50S 




$ 

FILE 

13,D3R,50S 




$ 

FILE 

14,D4R,50S 

a 




$ 

FILE 

15,D5R, 5iOS 




$ 

FILE 

H*,F6R,144R 




$ 

SYSOUT 

18 




$ 

SYSOUT 

19 




$ 

INCODE 

IBMF 




DATA 

21 

19. 




PARAM 

01 





CONSTR 

01 , 

4476. 

169. 



AREA 


21. 

45. 

270. 

315. 

$ 

ENDJOB 
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3.5 


PROGRAM OUTPUT 


The program output which would result from the setup 
described in Section 3,4.2 has been reproduced on the following 
pages. The printed output shows the area and the number of 
parameters used (Table 4), A tally of rev types and a list 
of revs for each type is printed (Table 5). The revs which 
were constrained are printed in Table 6, Tables 7 and 8 differ 
only when passes are edited because an insufficient number 
of crossings exists to solve for the requested set of parameters 
for each pass, E.g,, it is not possible to estimate both bias 
and slope parameters for a pass with only one crossing, and 
the pass is accordingly deleted. The allocation of memory 
for matrices and their starting addresses are shown in Table 
9, Table 10 shows the biases found (note that Rev 4476/Seg. 

169 was constrained) , Tables 11 and 12 show the corrections 
made at the crossing points (respective of the north- to- south 
in Table 11 and the south- to-north in Table 12) . 

Additional output has been included to show the effect 
that passes reaching maximum latitude have on the program 
i.e,, how the program handles passes that can be typed in 
more than one way (Tables 13-21), 

The data which is corrected and output on 9 -track tape, 
readable by the ECLIPSE S200, is described in Table 22. 
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SEAHT 


LEAsT'SQUaRES BETERMINATION of SkiOOThED sea surface HBIGHIS 

ca-oulated at intbrsecting points of crossing revolutions 


CAUIBRATICN area 
latitude 21,00 TO 45,00 

longitude 270,00 TO 315^00 

No'.oc parameters » i 

-NO, 8P SOUTR-TCi-NORTH REVS q 7 
-N0u.B,E-_Nt3RTH-T0«.50UTH REVS 


TABLE 4 



, ORIGINAL PAGE 
f OF POOR QUAII 


NORTh-TO-SOUTH 

REVOLUTIONS 

B6V NO 

SEt» NO 

4476 

169 

4334 

82 

4604 

249 

44^2 

160 

43^1 

117 

1974 

294 

1576 

462 

2102 

387 

1178 

104 

3552 

079 

1718 

81 

1988 

303 

south-to-north 

REVOLUTIONS 

8ev no 

SEC NO 

4610 

255 

4624 

278 

4482 

175 

4553 

218 

_2094 

382 

1710 

73 

21^1 

422 


TABLE 5 



»>»aofl#y5*«ttC0NSTRAINED REVS»********** 


-REV NQ. , SBG.NO. 
A476 ■ ■ 169 


45 . 

C^ 


TYPE PF constraint 
HEIGHT SLOPE CURVaTURE 


TABLE 6 



AflC CROSSIMO data INPUT FOR THIS R‘JN 




REVOwUTI 

JJOR’^HoSOOTh 

('TlV/SfcSJ 
4476 i69 

4476 J69 ' 

4476 j69 


4476 J69 ■ 

4334 89 

4334 Bg 

4334 02 

4334 02 

4334 02 

46q4 249 

46i)4 249 

49o4 249 

4Sj4 249 

46g4 249 

4$o4 249 

46n4 249 

4462 

4J’l 117 

1?74 294 

1?74 294 

1?74 2’4 

i?74 294 

i?74 294 

1?74 294 

lP74 294 

j976 4&2 

l976 462 

J.S76 462 

lS76 462 

i976 462 

ll76 462 

2102 S07 

2102 807 

2102 307 “ 

2102 337 

2102 307 

2102 807 

1170 t04 


OM NO 

south-north 
(REV/Seg) 
4610 255 

4624 270 

4402 i7| 

4553 210 

2j94 392 

1710 73 

4610 259 

4624 270 

4492 l75 

4353 218 

1710 73 

4610 259 

462 ^ 27* 

4192 i?9 

4953 2j.8 

2094 392 

l7i0 73 

fill 5?i 

4432 i78 

4610 259 

4624 270 

4402 i?5 

4553 210 

2o94 392 

1710 73 

2l5l 422 
4610 255 

4482 l75 

' 4553 21 * 

2C|94 332 

I'^io 73 

2151 422 

4610 259 

4482 179 

4553 218 

2094 332 

I7l0 73 

2151 422 

4610 259 



crossing points smoothed sea surface HkIGHT 

1,aTiTUd= LONGiTUpfc RORTh-SoUTH SOUTh-NOHTH 


(DEGRSES) (DEGREES! 
32,45963 29l. 04935 

35,09903 294,50023 

33,33520 292.52420 

34*63956 293,55636 

27,6723g 266,30274 

3o«50o33 29 d,446q6 

35,o7B0g 209,76633 

30,29404 292,52493 

35,09952 290.46118 

37,11724 291,49323 

33*31563 286.30353 

33,335oT 291.17394 

36,7oi4j 293,91276 

34,19322 291,84302 

35.4630'j’ 292.00094 

20,6437g 287,70/45 

31.49521 289,77(179 

26.14944 205,99549 

37.o9i9a 290.01010 

37,09222 206,75306 

33,3790a 29l,i3?6t 

36.74193 293,07820 

34,236i9 291,01462 

35.51037 292.04657 

28.69168 287.67375 

3l. 54072 209,73604 

26,l934t 285,96199 

34,23416 290.46656 

35,07230 291.14130 

36.31009 292.17323 

29,64004 267,00007 

32,43319 289.06383 

27,l92<»n 285,20925 

34,25205 290.45231 

39,00976 291,12/11 

36,33569 292,l59o7 

29,6602*' 286,?0630 

32,45202 289, 04^40 
27,21304 205,27465 

35.05774 209,80349 


(MeTERS) (METERS) 


-53. 

41369 

-53.21657 

«44 

48449 

-44,83744 

"6l 

29365 

-5l, 40305 

<'47* 

9820J 

-48.35399 

• 58, 

62893 

-57,23110 

• 58, 

16207 

•57,5/306 

• 5i 

44165 

-51,51245 

o46. 

23443 

-46. 0l5o6 

• 49, 

84351 

-49,84620 

-48,47016 

"46 

38004 

• 56. 

32589 

-56.58433 

• 52 


•53.0/022 

• 45 

•45,46606 

e§0 ■ 

79486 

-51.44217 

•48,47261 

s47 . 

71143 

-65 

68108 

-55,99393 

c56 

72089 

-57,42927 

i6n. 

69884 

-5i ,54139 

-47. 

02329 

-48,10359 

*•48, 

02431 

-48,74084 

• 52. 

756(51 

-53,ldo49 

c45 

1733.1 

-45.43630 

•50 • 

06500 

-5i, 42990 

• 46. 

20108 

■■ .40.53621 

-65. 

64155 

•56,03102 

• 56 . 

99460 

•57,40225 

•5i, 

02344 

-51,39341 

"Sg. 

99675 

-53*02664 

•5l I 

52079 

-51*52054 

c4B, 

26250 

-49,l97lt 

-55, 

10346 

-55 , 01010 

-56, 

91794 

-56.90405 

•5o, 

34713 

•50.47452 

.52, 

33057 

•52.90447 

• 5o . 

66233 

-51,40000 

-48, 

16755 

-49,26047 

•55, 

05340 

-55,05220 

■ 56 ) 

01059 

-56,96490 

• 5i* 

23919 

-50.52425 

• Si, 

47779 

-51,57277 


TIME FROM EPOCH 
NoRThaSOUTh SOUTh-NORTH 


(MINUTES) 
2144,82202 
9i4.9q9i6 
1732*10018 
1115*11140 
4378 , 14001 

3025,93799 

1062*11766 

340*31153 

1491,41762 

9o/. 62002 

27i7, 12609 
1783.27779 
182.41309 


1377,85480 

771.75292 

3979.45212 

2649,2633/ 


5133.9563/ 

535.40232 


532.98451 

1688,78401 


91,68554 


1284,28571 

679,60391 


3080.30493 

2552,05614 

5032*57422 

I202.3l4l9 

005.12060 

2ll,5o9o3 

3356.01572 

2051,03781 

4489,88655 


1478,91104 

1081*68381 

488,02697 

3632,77325 


2327,69366 

4766,73920 


1606,56949 


(MINUTES) 

2094.13519 

2053,33249 

2144.07568 

1714,92035 

2002. 44300 

197o,3o9o 1 
3329.28116 
3200.98703 
3357.28058 
2095,14792 
3249.20496 
25o5 ,9 oq 63 
2456.21045 
?548, 61484 
2100,57404 
2450.497/4 
2396,45407 
239i .24704 
3267.14432 
4311.45142 
2526.51672 
2455,56882 
2568,91663 
21?0. 40247 

2472,67120 

2417.68823 

2414,17776 

2929,42444 

2964.67441 

2513.41693 

2911,67218 

2034.34293 

2071.02502 

2937,67918 

2972,95511 

2521.43231 

2921,04380 

2043,66019 

2080,85709 

3319,29059 


TABLE 7 



1178 

^04 

4n24 

278 

30, 

27530 

292.54100 

1178 

104 

4402 

175 

35. 

873g7 

290.47023 

Il7'8 

‘ 104 

4553 

218 

37, 

09769 

291,51006 

1178 

S04 

" 2q94 

382 

30, 

55520 

206,33011 

;i78 

-104 

1710 

73 

33. 

29426 

200,40095 

1178 

104 

2l5l 

422 

20. 

l54ll 

284,62661 

3552 

879 

4610 

255 

33. 

39455 

291,12039 

3552 

879 

4624 

278 

36. 

75532 

293,66666 

3»&2 

879 

4482 

175 . _ 

34. 

25o5? 

291,00320 

3552 

679 

4553 

218 

35. 

52424 

292.03!!,i8 

3552 

879 

2C’'» 

332 

28. 

70793 

207,66233 

3552 

879 

l7l0 

73 

3i. 

55601 

269,72542 

3552 

879 

2l5l 

422 

26, 

21544 

285,95057 



4^10 

255 

3i. 

57794 

292,51490 


®1 

4< 2-4 

278 

S'?, 

o912? 

295,25341 

l?l8 


4482 

l75 

32. 

47i8(, 

293.1BV75 



4553 

218 

33. 

8o306 

294,22169 


®1 

2o94 

332 

26. 

69729 

289,04087 

l7l8 

®l 

I7l0 

73 

29, 

66o3l 

291,11102 

l7i8 


2l5l 

422 

24, 

Ul7p 

287,33712 

- = f 

8o3 

4610 

255 

29. 

71754 

293,877i9 

IV 86 

303 

4624 

276 

33. 

38137 

296.61583 

1?8B 

303 

4<82 

175 

5o. 

64354 

294,55204 

-p- lif® 

803 

4553 

21 » 

52. 

03658 

295,50401 

00 l ?88 

303 

2:94 

332 

24. 

64896 

290,41068 

i?88 

803 

I7i0 

73 

27, 

72268 

292,47371 

1?88 

303 

23 51 

422 

21, 

'97517 

268,69078 


-.46,46430 

=^0 . 023()4 

s«0.l4S>?7 
-34.91659 
-B5. 9^979 
oBl, 34346 
-*6.94230 
'"41,2021? 
s45 , l9866 
'>43.015o5 
-46,7550? 
-49,95370 
-*0 • 58ip3 
-52,79746 
o44 , 53096 
-5[| . 6476 q 
*47.08(157 
-58 I 2?546 
-57,16647 
*52,06040 
-52,42044 
-44 , q7 055 
-6fl , 8o 455 
■44.49924 
•58,42932 
’•55,49704 
•53,67839 


.46,03265 
-49.83742 
-48,50958 
-54,<’8i99 
•56,59882 
-51,75623 
-53.11638 
•45, 43241 
•5l , 430 0 7 
-48 , 553 03 
-56, g425o 
•57,49687 
•51,34337 
•52,65 i7i 
-45.o6i53 
•50.70387 
-48 , o48g6 
•58,05201 
-57,30741 
-52,86401 
-52, 65i5o 
-44,13325 
•5l , 7 o6 58 
-44.72337 

-58.6977S> 

-55,92331 

-54,06509 


67,42700 

12l6, 52719 
633,73615 
3723,01672 
2440.21732 
4838.60590 
1844.67471 
246,14857 
1439,02391 

834.60037 
4038.01172 
27q9, 48355 
5191,16931 
167?. 58382 
16 , 48501 
1253,44106 
6?6, 49030 
3939,24167 
2567,31482 
5127 ,530o9 

1959.16600 

243.44985 

1525,30772 

875,87438 

4300,31201 

2884,51242 

5523.58356 


3347,54004 
2885.77856 
3337,09067 
3238,70004 
3314,17050 
2533,40982 
2461,97546 
2575,68604 
2134,99152 
2480.18259 
2424.82361 
2421.995/0 
l6«l, 29201 
l669,4io60 
1738,45343 
1320,14997 
1553.52240 
l543, li702 

1458,96642 
814.55627 
861,44814 
886,32440 
49o. I8060 
614,65333 
647,5o87o 
485.65553 


TABLE 7 (Gont.) 



arc crossing Data USEO EOR ADJUSThlENTS 


revolution ho 
NORThoSOOTh SOI-TH-HORTh 


(REV/StG) 

(FEV/SEG) 

4476 

I69 

4610 

255 

4476 

169 

4624 

278 

4476 

l69 

4402 

175 

4#76 

169 

4653 

218 

4476 

169 

. 2C94 

302 

4476 

■ 169 

1710 

■' 73 

4334 

02 

4610 

255 

4S34 

82 

4624 

279 

4334 

82 

4402 

175 

4334 

82 

4553 

218 

4334 

02 

1710 

79 

46q4 

2<9 

4610 

255 

4604 

249 

. 4624 

278 

46q4 

249 

4402 

l75 

46q4 

249 

4553 

210 

46q4 

249 

2c99 

302 

46q4 

8<9 . 

I’lO 

73 

46o4 

249 

2l5l 

422 

4462 

160 

4553 

2l8 

4591 

117 

4402 

175 

•p* 

294 

4610 

255 

1?74 

2’4 

4624 

278 


2’4 . 

4482 

175 

1974 

394 

4553 

21» 

1974 

2’4 

2p94 

302 


a’4 

1710 

73 

l974 

294 

2151 

422 

1»76 

462 

4610 

255 

. 1>76. 

<62 

4402 

l75 

1.576 

462 

4553 

21* 

1576 

462 

2C94 

392 

1576 

462 

1710 

73 

l576 

462 

2151 

422 

2102 

307 

4«10 

255 

2102 

. 387 

4482 

175 

2102 

387 

4553 

2l8 

2102 

387 

2C94 

332 

2102 

397 

1710 

73 

2102 

307 

2l5l 

422 

1170 

104 

4610 

255 


CRoSSIMG points BHOOThED Sea surface heiqht time from epoch 

LaTiTUds longitude NORTh-SoUTH SOUTh-NORTh .. north-south S0 UTH“N0RT> 

IPEGREES) (DEGREES) (MeTERS) (METERS) (MINUTES) (MINUTES) 

32,45966 291,84935 o53,4l366 -53,21657 2144.62202 2094,13519 

35.899Q3 294,58823 -44,40449 -44,83744 5i4.9q9i6 2053.33249 

33*33520 292,52428 -»5i, 29365 -51,48305 l73?,l8oi8 2144.07568 

34,63356 293.55638 s47, 96283 -48,35899 1116,11148 l7l4.9?035 

27,67280 280,38274 _ -58,62896 -57,23110 4378,14001 2002,44380 

30.56033 29o,446o6 -50.i62()7 ' -57,57306 '’ 3025,9379? ' l97o,3o9oi 

35,d788c 289,78633 -5i, 44165 -5l, 51245 l882,ll766 3329,28116 

38,29404 292.52493 *46,23445 -46,015q6 340*31153 3200.98703 

35.89852 29o,46u8 •49,84351 -49,84628 I49l,4l762 3357,26050 


37,11724 291,49323 -48.38054 

33.31663 268,38353 -56,32589 

33.33507 291.17394 - 52 , 53 l 5 l 

36,70145 293.91276 •45.11979 

34,19322 291.64882 -5o,79486 

35,46365 292.B8b 94 »47,7ll43 

26,6437c 287.7 o/45 -55.68i86 

31,49321 289,77d79 •56.72oB9 

26,14844 285.99549 -6n, 69084 

37,69196 29o«8i8io -47,82329 

37.89222 203,75386 -40,02431 

33.37906 29l.l3?8l •52.756()i 

36,74193 293,87028 -45.17331 

34,23519 29i, 81462 .So.aBSRe 

35,51037 292.04657 •48,20100 

28,69168 287,67375 *55,64155 

3i>34072 269,73664 *56,99466 

26,19341 265.96199 *5i, 02344 

34.23416 290-46656 •52,99675 

35 ,o 723 c 29 i , 14130 « 5 i, 52 o 79 

36.31089 292,1/323 ■48.26256 

29,64004 287.00087 -55,lo346 

32,43319 289,06383 •66,9i784 

27,19248 285.20?25 *5o.347l3 

34,25205 290.45231 -02,33857 

.35*08976 291,12711 *50,66233 . 

36,33569 292,15?o7 -48,16755 

29,66026 286,?8638 -55,0534(J 

32,45202 289,04948 •56,6t059 

27,21384 265.27465 .“5i,239i9 

35.05774 209,00349 -5i, 47779 


-40,47816 9o 7. 62082 2095.14792 
-56,58433 27i7. 12689 3249,2o496 
-53 ,o 70Z2 1783,27779 Z5o5,9oo63 
-45.46606 182.41369 2436,21045 
-51.44217 1377,85480 2548,6i484 
-48.47261 771,75292 2103,57404 
-55.99393 3979,45212 2450.497/4 
-57,42927 2649,2633/ 2396,45407 
-51.54139 5133*9q637 239i,247o4 
-48,18358 535,48232 3267.14432 
-48,74084 532,90451 4311,45142 
-53*10049 1686,78401 2526,5i672 
-45,43650 91,68554 2455,56882 
-5l, 42990 1284,26571 2568,9i663 
-48,53621 679,60391 21?8,4o247 
-56,03182 3880*30493 2472,67120 
-57,46225 2552.85614 2417,68823 
-51.39341 5032*57422 24l4. 17776 
-53.02664 12o2*3l4l9 2929.42444 
-51*52854 8o5. 12060 2964,67441 
-49,1?7 ii 2i1.50903 25i3.4i 693 
-55,01818 3356,01572 29ll. 67218 
-56,90485 2051,03781 2034,04293 
-50*47452 4489,88635 2071,02502 
-52,90447 1470,91104 2937,87910 
-61,48080 1001,68301 2972,95511 
-49,26047 488,02897 2521,43231 
-55,05228 3632,77325 2921,04388 
-56.96498 2327,69366 2843.66019 
•50*52425 4766,73920 288o,857o9 
-51,57277 1606.56949 3319.29059 
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1178 

J04 

4624 

278 

1178 

104 

4432 

175 

1178 

104 

4553 

21« 

1176 

104 

2q94 

302 

1175 

104 

1710 

73 

1178 

104 

2151 

422 

3552 

679 

4610 

25S 

3552 

879 

4624 

278 

3352 

8 79 

4482 

175 

3952 

879 

4553 

21* 

3952 

879 

2094 

302 

3952 

879 

1710 

73 

3952 

679 

2151 

422 

ljl3 

Si 

4610 

255 

l7lS 

-Si 

4624 


l’l6 

" - 8.1 

44(8^2 

175 

l^lS 

'8~i 

4553 

21 “ 


81 

2q94 

302 

I’lS 

61 

I7l0 

73 

thf 

81 

2191 

422 


503 

4610 

259 

1?08 

303 

4624 

278 

1?8B 

303 

4482 

175 

1?86 

303 

4553 

2l» 

l?b8 

303 

2c94 

392 

1?B8 

303 

1710 

73 

lf88 

303 

2151 

422 


38,27536 292.>'>158 
35,B73q7 29o, 47(J23 
37,09769 29i, 51006 
30.55623 286,33?H 
33,29428 288.40095 
2P.15411 264,62661 
33.39451 291,1283’ 
36,75532 293,86686 
34,25052 29l. 80320 
3 , 5,52424 292.83518 
2fl*, 70793 287,66233 
3li556oi 269,72542 
26,21544 285.9t)057 
31.57794 292.51498 
'•35,0’123 295,25341 
32,47i8o 293,i8?75 
33,Bo3o6 294.22169 
26,69729 269,04887 
29,66031 291.11182 
24.11172 287.33712 
29.7i754 293.877i9 
33,38137 296,61583 
30,64554 294,55204 
32,03658 295,58401 
24,64396 290,41068 
27;72268 292,47371 
2l,975i7 288, 69^76 


.46,46430 

-46, 

"5o. 02304 

-49. 

•48,14527 

• 48, 

"54 , 9i659 

-54 , 

=55,9i979 

-56, 

•51,34346 

-5l. 

•,46,94230 

•53. 

'’41,20212 

-45, 

•45,19866 

«5i , 

•’43,01505 

■ 48 , 

•46 , 755o7 

• 56, 

•49,95378 

•57. 

’’40,58103 

•5l , 

-52,79746 

-52 • 

"44. 53096 

-45 , 

•5o , 64760 

•5o , 

-47,86057 

-48 , 

o58 , 2?546 

-58, 

-57,16647 

-57, 

-52,86048 

•52, 

-52.42044 

•52, 

•<14 , 07q55 

-44, 

-5q, B o455 

•5l, 

-44,49924 

-44, 

"58,42932 

-58, 

-55,49784 

• 55 , 

- 53,67839 

-54, 


03265 67,42700 
33742 1216,52719 
50958 633,73615 
93i99 3723,01672 
59882 2440,21732 
75623 4836,60590 
11638 1844,67471 
43241 246,14857 
4300? 1439.62391 

55303 634.60037 
04250 4038,01172 
49687 27o 9, 48355 
34337 5191,16931 
65i7l 1672.58382 
06i53 16 . 48501 
70387 1253.44106 
048q6 626,49030 
0>201 3939.24167 
30741 2567,31482 
86401 5127,53009 
65i5o 1959,16600 

13323 243,44985 
70658 1525,30772 
72337 375,87438 
69779 4300.31201 
92331 2884,51242 
06509 5523,58356 


3191,98083 
3347.54004 
2885. 77856 
3337,09067 
3236,70004 
3314.17050 
2533,40982 
2461,97546 
2575.68604 
2134,99152 
2480,18259 
2424.82361 
2421,995/0 
16Si,292U1 
l669, 41060 
1738.^15343 
1320.18997 
1553,52240 
l543.1i702 
1450,96642 
8lil, 55627 
661,44814 
886,32440 
490,1B060 
614.65333 
647,5o670 
485,65553 
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STIRTING AD0B6SSE5 OF ARRAYS 

SHREy 1267CD0GOflf>0 
WSHty i2«71f003000 
SUEl 12472^030000 
SMH2 12703«009000 
fiEUTA 127054030000 
P81AS 12707AOOOOOO 
N?TR J27104000000 
8R?AY 127106000(T08 
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TABLE 9 



REV NO 
44!t6 
44^4 
46p4 
44^2 
43»l 
19?4 
1576 
2102 
ii?a 

3552 

1718 

1990 


SSG NO 
169 
82 
249 
160 

- 11 ? 
294 
462. 
387 
104 
B79 
01 
303 


NORTH-TO^SOUTy RCvS 


COEFncIENTS 

70i?l 66433; 
0,2fi0 451'5 = 
0.7391667= 
0,43 81491 = 
0,8481592= 
0,5389130E 
0.2523978= 
0 .4ii47830 = 
0.2H8i538= 
0 > 7;'54924 = 
0.2282054= 
0.5026312= 


14 

00 

00 

00 

CO 

00 

00 

DO 

00 

01 

DO 

DU 


REV NQ 

SOUTH- 

tojinqrth revs 

SSG NO 

CSEFT1CIENT5 

4610 

255 ‘ 

0.2682993= CC 

4634 

278 

0.5*147675= 00 

4482 

175 

0. 1236266= CO 

4553 _ 

218 

0.5/849475-01 

20?4 

332 

0.1714784= 00 

i7i;o 

73 

0.5947062=-01 

2151 

422 

“0.2476388= 00 


Ln 

IsJ . 


TABLE 10 



ReVOl.UT10NS 
N« S/N 

CR0SSIN3 PTS; 

latitude longitude 

{DEGREES} (DEG5EES) 

BUSES 

N/S N/S 

JIETERS} (METERS! 

SEA SuftpA^^^BilGHTS 
S/N N/S 

liMETERS) (METERS) 

SEa SU^pAci^iilGHTS AVERAGE 
S/N N/S (ADJUSTED) 

(METERS) (METERS) (METERS) 

RESIDUALS 

N/S 

.(METERS) 


4615 

32,459675 

291,849334 

<•0,0000 

0<2I383_‘ 

_ "53,4137 

-53.2166 

-53,4137 -53.4849 

- 53.4493 

0,0356 

4874 

' 4624 

85,399035 2*4,563226 

• 5 U 000 

0»5fi46; 

■ 1 .- 44,4345 

-44,8374 

744.4845 -45,4222 

-44,9533 

0.4639 

_#S7^ . 

...4402 _ 

„33, 335204 

2?2,524231.. 

.-0,0000. 

.1_0,1?36: 

Hr5l,2937 

•51,4830 

. -51,2937 -51,6067 

-51,4502 

0*1565 

487$ 

45SJ 

34,630364 

213,556377 

-9. OQOO 

0.DS78: 

:7-*7,9328 

-48,3590 

•47,9828 -48,4168 

-48,1996 

0.2170 

4874 

20’4 

27,672798 

268,532736 

-oioooo 

” 0.1715' 

-58.6290 

-57,2311 

-58,6290 -57;4o26 

-58,0158 

-0,6132 

. 4S7$ 

.1710 _ 

3g, 330333 

2*0,446030 . 

,-0.0000. 

0 . 0595 . 

„ .-58. 1621 

-57,5731 

-58,1621 -57,6325 

. -57,89/3 

-0.2646 

4434 

4610 

35,j788o4 269,736335 

0.2009 

0.2B83 

-51,4417 

■51,5124 

-5l,7o25 -51.7807 

-51,7416 

0, 0391 

4534 

4624 

38,294036 

2?2. 524925 

0,26(j9 

0i5?48 

1. -46,2345 

-46,0191 

-46,4953 -46,5998 

-46,5476 

0, 0523 

4534 

440J_. 

35,898524 

2*0.461178 

. 0.2009, 

0'1?36„ 

" -49,8435 

•>49,8463 

-5q,lo44 _ -49,9699 

-50.0371 

-0,0672 

4534 

456J 

37,117242 291,493225 

O,20Q9 

0.0578 

-48, 38o5 

-48.4782 

-48,6414 -48.6360 

-48,5867 

-0.0527 


1710 

33.316635 

268.883526 

0.2009 

0.0595 

-56,3259 

-56,5843 

-56,5867 -56.6438 

-56,6153 

0,0285 

^ 44o4 

.*0*P 

33,336066 

2fl, 873943 

0,7392 

0.2883. 

] -52,5315 

-53,0702 

753,2707 -53,3365 

-53,3046 

_ 0. 0339 

4404 

4624 

36,7q1446 

2t3, 912762 

0,7392 

0,5648 

-45,1198 

-45,4661 

- 45,8590 -46,0508 

- 45.9549 

0.0959 

44o4 

4462 

34,193219 2 * 1.648816 

0,7392 

0.1336 

" -50.7949 

-51,4422 

“51,5341 -51,5658 

-51,5499 

0.0159 

._44Q4 

„*5i»S .. 

35.468647 

2*2.880936 

0,7392 

0.0578_ 

-47, 7114 

_ -48 , 4726 

-48,4506 -48.5305 

-48,4905 

0.0399 

4004 

2094 

28.643703 

267,707451 

0,7392 

0.1715' 

■ ■ >55. 6819 

-55,9939 

-56,4211 -56.1654 

-56,2932 

-0.1278 

4504 

1710 ■ 

31,495209 

269,770786 

0,7392 

" 0.0595 

' „-55.7209 

-57,4293 

-57,4601 -57;4887 

-57.4744 

0,0143 

,.440* 

...2l5l_ 

„ 26,148443 

265,f95491„ 

__0.7392. 

(,0,2876 

-50.6968 

_-5l,54t4 

-_-51,43ao„ -51.2937„ 

.. -51,3659 

.-0,0721 

. 4468 

.4553 

37,891960 

2*0.816096 

0.4181 

0,0578. 

;■ -47,6233 

-48,1836 

-48,2414 -48,2414 

-48,2414 

0, 

■ 439* 

■ 4402 I 

37^892282 

268,153857: 

".“o'|8482 

r"0tl336; 

”■.3*46,0243 

-48,7488 

-48,8725 -48,8725 

-48,8725 

0, 



TABLE 11 



J?74 

4610 

33,379863 

791.139809 - 

0.5389 

0.2^83 

-52*7560 

-53,1005 

-53,2949 

-53.3668 

“53.3319 

0.0369 

3L974 

4624 

36,741929 

293,978277 

0,5389 

. 0.5848 

-45,1733 

-45,4363 

-45, 7^22 

”46,0211 

- 45,8666 

0 • 1544 

1?74 

4482 

34. 236194 

29i,8146i7 

0.5389 

0.1236 

-50,6659 

-51,4299 

”51, 4c48 

”51.5535 

-51,4792 

0>0744 

1974 

4553 

35.510322 

292,046565 

0,5389 

0,0578 

-48.2010 

-48,5362 

•48.7?99 

-48,5941 

-46.6670 

-0«0729 

. 1974 

20’4 

28.691679 

287.673762 

0,5389 

„ OiUlS. 

-55, Ml6 

-56,0318 

-56.1F05 

-56.2033 

-56.1919 

0.0114 

1974 

1710 

..31.540717 

299,T36P35 

0,,53B9 

. 0*0595 

-56.9947 

-57,4823 

-57,5336 

-57.5417 

-57,5376 

0.0041 

1974 

2l5I 

26.198408 

285.961994 

0,5389 

50-2P76 

-51. 0234 

-51.3934 

-51,5624 

-51.1458 

-51,3541 

-0.2083 

157« 

4610 

34.234157 

290.466*64 

0,2524 

, 4 0.2S83 

-52.9967 

-53,0266 

-53.2*i9l 

-53,2949 

-53,2720 

0-0229 

15.7a 

4482 

35„C72296 

291,^.415.00 

0,2524 

0,1?36 

“5l , 5208 

-51,5265 

-51,7732 

-51,6522 

-51,7127 

- 0 " 06o5 

1576 

4553 

36.318888 

2’2.173225 

0,2524 

0*0576 

-48,2626 

-49,1971 

-48,5i50 

-49.2550 

-48,8850 

0-3700 

1576 

2g94 

,, 29.64004c 

287, (?(jo8a6 

0,2524 

. ,0*1?15 

-55 , ] 035 

-55, 01»2 

-55,3559 

-55:1097 

-55,2728 

”0.0831 

1576 

1713 

32.433169 

2?9,,463P31 

0.2524 

.0*0595 

-56,9178 

.-56,9848 

-57,1702 

-57 , 0443 

-57.1073 

”0-0630 

1576 

2l5j 

27.192479 

285.,28 9?4 6 

0,2524 

JO* 2 fi 76 

-50,3471 

-5o, 4745 

-5o,5?95 

-50,2269 

-50,4132 

-0. i663 

21Q2 

4610 

34.^252049 

2?0, 452512 

0,404B_ 

0>2®83 

-52.3386 

-52.9045 

-52,7434 

-53.2528 

-52,9981 

0-2547 

2102 

4482 

35,089757 

29i,f27jo6 

0.4046 

0*1236 

-50 . 6623 

«6l , 48o8 

-51,0671 

-5J , 6c44 

-51,3358 

0-2687 

2102 

455>3 

36.3.35685 

2?2,159('65 

0,4048 

0*0578 

-48 . J676 

-49 ,26o5 

-48,5723 

-49,3183 

-48,9453 

0.3730 

2102 

2(J94 

29.660255 

?86, 986:182 

0,4043 

_.0'1?15 

“55,0534 

-55,0523 

-55.4582 

-55,2238 

-55,3410 

”0-1172 

21CH 

1710 

32.452020 

289, 849.*80 

0,4043 

0,0595 

-56,8106 

-56,9650 

-57,2i54 

-57,0245 

-57,1199 

-0.0955 

2102 

, 2li>I 

27.213644 

285 ,, 274654 

0,4048 

**0*2?»76 

-51.2392 

•5c. 5242 

-51, 6<i40 

-50.2766 

-5 0 , 96 u3 

-0-6837 

i;76 

4610 

35.057742 

289,9q3«90 

0,2882 

0*2<5B3 

-51,4778 

-51,5728 

-51,7659 

-51.6411 

-51,8035 

0-0376 

1179 

4624 

38.275364 

292,541580 

0.2882 

0*5848 

_ -46.4643 

-46,0327 

-46,7525 

-46,6174 

-46,6849 

-0-0675 

1?7§ 

j4482 

^ 35.878066 

290,478233 

0,2862 

0,1236 

j50. 0230 

-49,8374 

-50,3il2 

-49.9610 

-50,1361 

-0-1751 
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iX76 

4553 

37.097687 29i,5ioQfi3 

0,2882 

- 0'0S78 

.-48,1493 

•48,5p96 

-48,4334 

-46,5674 

-48,5004 

0.0670 

_U7^ 

20’4 

_ 30.556260 286,?38j.o8 

0,2882 

0.1Z15 ■ 

■ -54, “166 

-54,9820 

-55,2047 

-55,1535 

-55.1791 

-0.0256 

i?79 

1710 

33.294276 298,0OO'M6 

0,2082 ■■ 

■ 0,0595 ■" 

-55,«198 

-56.5988 

•56,2o79 

-56,6583 

-56,4331 

0.2252 

1^70 

2l5i 

28.154114 284.626.So6 

0,2082 

40*2j!76 

-51.3435 

-51,7562 

-51,6316 

-51,5086 

-5l,57oi 

-0,o6i5 

, J558 

_ 4610 

33. 394506 29l.j2BS?l 

_7,2549_ 

_„0i29B3__ 

..-46,9423 

-53,1164 

-54,1972 

-53,3847 

-53, 7909 

-.0.4063 

?S52 

4624 

36.755324 293,866856 

7,2549 ■ 

0,5640 

-41.2021 

-45.4324 

•48,4570 

-46,0172 

-47.2371 

-1,2199 

3552 

4482 

34,250523 29i,8o33o4 

7.2549 

0*1?36’; 

“-45,1987 

-5i,43o1 

-52,4536 

-51 ,'5537 

-52,0036 

-0.4499 

3552 

4553 

35,524237 292,635178 

.... 7,2549__ 

0*0598„ 

.-43,3151 

-4B.5530 

.-50,2700 

-48.6109 

-49,4404 

_,.-0.8295 

3552 

2|}74 

20,7o793o 287,662331 

7,2549 

0,lZl5 

-46,7551 

•56,0425 

-54,0100 

-56,2140 

-55,1120 

1.1020 

3552 

1710 

31,556008 289,)[25SiB 

7,2549 ■ 

0.0595 

-49,9538 

•57.4969 

-57.2007 

-67-. 5563 

-57,3825 

0.1730 

3552 

2151 

26.215442 285,850865 _ 

. 7,2549 _ 

".J0-2I»76_ 

_-40.SBl0 

-51,3434 

-47 1 6360 

-51.0957 

-49,4658 

1,6299 


4610 

31.577943 292,514976 

0,22a2~ 

0,2883 

-52,7975 

-52.6517 

-53 , 0257 

-52,9200 

-52,9720 

-0.0528 

I7i9 

4624 

35.091232 295,jS31o7 

0,2282 

0,5§48 

-44 , 3310 

•45, 06i5 

-44,7592 

-45,6463 

-45,2027 

0.4436 

iZl9 

4482 

32,471801 293.189754 

0,2282 

_ 0'1?36 ■ 

-50-6476 

-50 , 7039 

-50,8758 

-50.8275 

-5o,85l7 

. -0.0242 

I7l8 

4553 

33,903059 ?94,g2l487 

0,2282 

0,0578 

-47,5806 

-48,0401 

-48,lo08 

-48,1059 

-48,1073 

-0. 0014 

I7l9 

20’4 

26,697291 289,J48970 

0.22B2" 

0.1Z15 

-58.2255 

-58,0520 

758,4537 

-58,2235 

-58,3306 

-0.1151 

1710 

. 1710 

29.660306 29l,fn8z4 

. 0,2282,„ 

.„.0i0595 „ 

_-57,l665_ 

-57,3074 

. -57,3947 

-57,3669 

-57,3000 

... -0.0139 

171? 

21&1 

24,111718 287,337124 

0,2282 

|*0,2t!76 

-52,8605 

-52,8640 

•53,0387 

-52.6164 

-52,8525 

-0.2362 

i?a? 

4610 

29,717535 2?3,6j7193 

0,5026 

0,2883. 

-52,4204 

-52.6515 

•52,9231 

-52.9198 

-52,9214 

-0.0016 

_ 1 ? 8 | 

, 4621( 

„ 33.381367 2?6.6i5329 

0.5P26.J;; 

I_0!5648'_ 

..-44,0706 

•44,1333 

-44,5732 

-44,7l8o 

-44.6456 

0.0724 

,l?a? 

4482 

30.648543 2?4, 552044 

0,5026 

0,1?36 ■ 

-50 , 8046 

•51,7066 

-51,3072 

-51,8302 

-51,5687 

0.2615 

1980 

4553 

32,036575 2?5,584(Jll 

0,5026 

1 0.0578 

■-44.4992 

•44.7234 

-45,0019 

-44,7812 

-44,8915 

-0.1103 
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20^4 

24.648965 

2*0 •'•10675 

0,5026 

0.1U5 

-58.4293 

-58,6978 

-58,9320 

-58.8693 

-58,9006 

-0.0313 

. 1?09 

1710 

27,722684 

2*2. <73713 

0.5026 . 

0.0!|95_ 

...-55,4978 

-55.9233 

"56,0005 

-55.9828 

-55,9916 

-0.0068 

i?ae 

2l5J 

21.975166 

208,6*8776 

0,5026 

10.2976' 

_-53.6784 

-54,0651 

-54,1810 

-53.8175 

-53.9992 

-0.l8l9 

MBAN OF 

RESIDUALS 

iQ. 09000 










STftN?ARD„ DEtflATJON 0,97309 


tn 

O' 


TABLE 11 (Cont.) 



A PRIORI - ADJUSTED 

REVQIUTIONS crossing pts, buses sea surface heights sea surpace heights average residuals 



N/S 

latitude 

(DEGREES! 

LQHGITUDE 

(0EG8EES) 

S/N 

JvtETERS) 

S/N 

(METERS! 

N/S 

(METERS) 

S/N 

(METERS) 

N/S 

(METERS) 

S/N 

(MEIERS) 

ADJUSTED) 

(meters) 

S/N 

(meters) 

^fiie 

4476 

32.459675 

2?ii$49j54 

0;2683 

(to 

0800 

”53 , 4166 

•53,4137 

•53,4849 

-53.4137 

- 53.4493 

-0-0356 


4334 

35.078004 

2«9,J86335 

0,2603 

0 

280* 

"5l, 5l24 

-51,4417 

-5i, 78 q7 

-51.7025 

-5l, 7416 

-0*0391 

46lfl 

_46U4_ 

33,336068 

291,173943 

o;2683 

0 

T'592„ 

_ -63,0702 

•52,5315 

-53,3385 

-53.2707 

_ -53.3046 

-0*0339 

4fllS 

1974 

33,379663 

291.J393o9 

O', 2683 

0 

5389 

-53.1005 

/ 

-52.7560 

-53,3608 

-53,2949 

-53,3319 

-0*0369 


1576 

34,234157 

29q, 166564 

0.2683 

0 

2524 

-53,0266 

-52,9967 

-53,2949 

-53,2491 

-53,2720 

-0*0229 

4^18 

, 2102 

,.34,252049 

2?0. 452312 

0.2683 

„ 0 

4p48__ 

_-52,9S45 

•52,3366 

•53.2520 

.-52,7434 

-52.9981 

-0 * 2547 

46ie 

1178 

35,057742 

209.803490 

0',26B3 

0 

2S62 

-51,5728 

-51,4770 

-51.0411 

-51,7659 

-51,6035 

-0*0376 

4i»ie 

3552 

33,394506 

291.128591 

0.2683 

■ 7 

2549 

-53,1164 

-46,9423 

-53,3847 

-54;l97a 

-53.7909 

0*4063 

4«10 

1716_ 

31,577943 

2»2. 514976 

. 0‘.2683 

0 

2?82_ 

_-52,65l7 

-52.7975 

.-52,9200 

-53,0257 

_. -52.9728 

0*0520 

46ie 

1988 

25,717535 

293,877193 

0'.2683 

0 

5p26 

-52,6515 

-52, 4204 

-52,9i98 

-52,9231 

-52,9214 

0. 0016 

4624 

4476 

35.699035 

294,588226 

0,5048 

40 

opoo 

-44,0374 

-44,4045 

•45,4222 

-44; 4845 

- 44.9533 

-0*4689 

tn 

•<1 462* 

4334 

36,294036 

2?2,5249z5 

0,5848 

_._o 

2flo9_ 

-46,0151 

-46,2345 

■46 , 5990 

-46,4953 

-46,5476 

-0*0523 

4624 

46u4 

36,701446 

293.912762 

0,5846 

0 

7?92 

-45, 4661 

-45,1198 

-46, 0500 

-45.8590 

- 45,9549 

-0* 0959 

4624 

1974 

36,741929 

293.676277 

0;5B48 

0 

5589 

-45.4363 

-45,1733 

-46,0211 

-45,7122 

-45.0666 

-0 . 1544 

4624 

1170 

, 36,275364 

292.541580 

. 0,5048 

0 

2802. 

_-46,q327 

-46,4643 

-46,6i74 

-46.7525 

. -46,6849 

0* 0675 

4624 

3552 

36.755324 

293.866356 

O'.504B 

7 

2549 

-45,4324 

•41,2021 

-46,0172 

-48:;457o 

-47,2371 

1* 2199 

4624 

1718 

35,091232 

295,853407 

0,5040 

0 

2?82 

-45,0615 

-44,5310 

-45,6463 

-44,7592 

-45,2027 

-0.4436 

4624 

_i«J88 

33,381367 

296,615929 

O', 5043 

0 

5 e 26.1 

H^44.1333 

.-44,0706 

-44, 7i80 

-44,5732 

-44,6456 

-0.0724 

4462 

4476 

33,335204 

292,524281 

0.1236 

90 

opoo 

-51.4830 

-51,2937 

-51 , 6q67 

-51,2937 

-51,4502 

-0*1565 

_46a2 

4334 

35,898524 

2*0.861178 

Otl236 

0 

2809 

.-49,6463 

-49,8435 

-49,9699 

-50,1044 

-50,0371 

0.0672 
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4604 

34.193219 

?9l. $46616 

0.1236 

0*7592 

"51.4422 -50,7949 

-51,6658 

-51.5341 

-51,5499 

-0-0159 

J«8H 

4391 

.37,692222 

298.753357 

.0,1236 . 

._0< 6S82 _ 

.-48,7408 -48,0243 

-48,8725 

-48.8725 

-48,8725 

0. 

MSB 

1974 

34,236194 

29i,8i46i7 

0,1236 

0,5§89 

-51.4299 -50,8659 

-51,5535 

-51.4048 

-51,4792 

-0 . o744 


V57 6 

35, 072296 

291,141300 

0,1236 

0*2524 

-51,5285 -51,5208 

-51,6522 

-5177732 

-51.7127 

0. 0605 


2102 

35,089757 

2?1, 127106 

0,1236 

0>^P4B 

-51,4808 -50.6623 

-51, 6q44 

-51.0671 

-51,3368 

-0.2687 

<*162 

1178 

35.873066 

290,478233 

0 1 1236 

0*2682 

-49,3374 -50.0230 

-49,9610 

-50,3112 

-50.1361 

0.1751 

4fl08 

3552 

34 .250523 

29i,$o72o4 

0.1236 

7.2549 

-51,4301 -45,1987 

-51.5537 

-52.4536 

-52,0036 

0.4499 

4562 

,1718 

32.471801 

2?3, 169754 

0,1236 

0*2?82 . 

-50,7039 -50,6476 

-50.8275 

-50,6758 

-50.8517 

. 0.0242 

4582 

1988 

3q, 649543 

1 

294,552044 

0,1236 

0,5026 

“5l,7066 -50,8046 

-51,8302 

-51.3072 

-51,5687 

-0.2615 

4553 

4475 

34,638564 

293,556377 

0.0578 

«0< 0000 

-48,359(1 -47,9828 

-48 , 4l68 

-47,9828 

-48,1998 

-0.2170 

4U53 

4334 

37.117242 

^9i,49r2v5 

0,Q578 

0'26o9 

-48.4782 -40,38o5 

-48.5360 

-48,6414 

-48,6887 

0.0527 

4^53 

46i<4 

35.468647 

?92,68p936 

0.Q578 

0*7592 

-48,4726 -47,7114 

-48,6305 

-48,4506 

-48,4905 

-0.0399 

4553 

4462 

37.891960 

29o,$lP0’6 

0,0578 

0,4181 

-48.1836 -47,8233 

-48,2414 

-48.2414 

-48,2414 

0. 

4553 

1974 

35.510372 

292,64C5ft5 

0.0578 

0,5589 

“46.5362 -48,2010 

-46 , 6941 

-48,7399 

-48,6670 

0.0729 

4553 

1576 

3/.. 318888 

292,173225 

0.0578 

0.2524 

-49.1971 -48,2626 

-49,2550 

-40 , 6l5 0 

-48,8850 

-0.3700 

4?63 

2102 

56,335685 

292.159065 

o;0378 

0,4048 

-49,2605 -48,1676 

-49,3i83 

-48,6723 

-43.9453 

-0,3730 

- 4553 

1178 

37,097687 

??l>5tP063 

0i0578 

0*2082 

-48.5096 -48.1453 

-48,6674 

-48,4334 

-48.6004 

-0,0670 

4553 

3552 

35,524237 

292,835175 

0i0578 

7,2549 

-48,5530 -43,0151 

-48,6109 

-50.2700 

-49,4404 

0.8295 

4§53 

1713 

33,803059 

294,221687 

O', 0578 

0,2?82 

-48.0481 -47,8806 

■48, 1059 

-48,1088 

-48,1073 

0.0014 

__4553 

. 1988 

32.038575 

295,'?34oU 

0,0578 

fl.5p26_ 

_-44,7234 -44,4992 

-44,7812 

-46,0019 

-44,8915 

0.1103 

ap94 

4476 

27.672798 

298,882736 

0;i7l5 

40,0000 

-57,2311 -58,6290 

-57 , 4o26 

-53.6?90 

-58.0168 

0.6132 


4604 

23,643703 

287,707491 

0,l7l5 

0.7592 

-55.9939 -55.6819 

-56 , 1654 

-56.4211 

-56,2932 

0.1278 
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ap94 

1974 

2«, 691679 

287,673792 

0il7l5 

0.5389 

-56.0318 

-55,6416 

•56,2033 

-56,1805 

-55.1919 

-0.0114 

2Q94 

1576 

29,640040 

287,500866. 

0il7l5 

-0.2524. 

-55.0182 

-55.1035 

_-55,1B97 

-55.3559 

-55,2728 

0.0831 

2594 

2108 

29,660255 

286,936382 

0,17i5 

0.4348 

-55.0523 

•55,0534 

-55,2238 

-55,4582 

-55,3410 

0.1172 

2P94 

U70 

30.556280 

286,?35io8 

01i7i5 

.0.2802 

-54.9020 

-54.9166 

-55,1535 

-55.2047 

•55,1791 

0.0256 

..ap’t.. 

3552 

23.7o793q 

287,662331 

_0’tl7i5. 

7., 25 49 

_-56,0425 

-46.7551 

. -56,2140 

-54i0100 

•55,1120 

-1.1020 

209<J 

1718 

26,697291 

289,o4'4870 

0.1715 

0.2282 

-56,0520 

-58.2255 

-58,2235 

-58,4537 

-58,3386 

0.1151 

299| 

1968 

24.646965 

290.410675 

0.17i5 

0.5826 

-56,6978 

-58,4293 

-50,8693 

-58.9320 

-58,9006 

0.03i3 


4476 

_ 30.580333 

2*0.946q60 

..010595 

•0.0P00_.. 

.-57.5731 

. -58,1621 

-57 ,6325 

. -56,1621 

-57,8973 

0.2648 

l7i0 

4334 

33.316635 

298,383526 

0.0595 

0.Zfi09 

-56.5843 

•56,3259 

-56,6438 

-56,5867 

-56,6153 

-0.0285 

1710 

45u4 

31.495209 

289,773786 

01o595 

‘ 0.7392 

-57.4293 

-56.7209 

- 57 , 4887 

-57,4601 

-57,4744 

-0.0143 

.l7l0 

1974 

_ 31.540717 

289,736835 

0.0595 

0.5389__ 

-57,4823 

-56,9947 

-57.5417 

-57:5336 

-57,5376 

-0,0041 

I7i0 

1576 

32.433189 

2»9, 46,3831 

010595 

0.2524 

-56,9048 

-56,9178 

“57 . 0443 

-57,1702 

-57,1073 

0.0630 

on 1713 

to 

2102 

32.452020 

289.943400 

0.0595 

0.4848 

-56,9650 

-56,6106 

-57,0245 

-57.2154 

-57,1199 

0.0955 

„ I7l8 

1178 

33.294276 

283,40)948 

, 0,0595 

.0.Z882__ 

-56.5988 

-55,9198 

-56,6583 

-56.2099 

-56,4331_ 

-0.2252 

,17ih 

3552 

3t.556ofi8 

289,7254x8 

0. 0595 

7,2549 

-57.4969 

-49.9538 

-57,5563 

-57,2087 

-57,3525 

-0.1738 

I7l3 

1718 

27.660306 

291.111824 

0,0595 

0.2?82' 

-57.3074 

-57,1665 

-57,3669 

-57.3Q47 

-57,38o8 

0,0139 

..1710 

1988 

27,722684 

292.473713 

0.0595 

O' 5026. _ , 

-55,9233 

-55.4978 

-55,9628 

-56,0005 

-55.9916 

0.0088 . 

a?5i 

4604 

26,148443 

?85, 995491 

-0.2476 

0.7392 

-51.5414 

-50.6988 

-51,2937 

-51.4360 

-51,3659 

0.0721 

2151 

1974 

26,198408 

285,961994 

-0.2476 

0,5369 

-51-3934 

-51,0234 

-51,1458 

-51.5624 

-51,3541 

0.2083 

_2tSl 

1576 

. .27,192479 

?»5, 299246 

-0.2476 

. .0* 2?24 

-50.4743 

-50 , 3471 

-50,2269 

-50,5995 

-50,4132 

0.1863 

a;3i 

2102 

27.213844 

285,274654 

-0.2476 

0.4848 

-50.5242 

-51,2392 

-50,2766 

-51.6440 

-5o,96U3 

0.6837 

3151 

ll78 

28.154114 

284,626606 

-0,2476 

0.2882' 

-51,7562 

-51.3435 

-Sl,5o86 

-51,6316 

-51,5701 

0.0615 


ORIGINAL PAGE IS 
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3;si 

3552 

26,215442 

?85,?5o565 

-0.2476 

7.2549 

-51,3434 

•4o.58ig 

-51,0957 

-47,8360 

“49,4658 

-1.6299 


__1718 ^ 

2-<, 111718 

287,337124 _ 

■<0!2476 

* 

0.2?82 

-52. 8640_ 

-52,8605 

-52, 6i64 

-53, 0887 

-52,8525 

0.2362 

tX5i 

1998 

21.975166 

288,6?3776 

"0,2476 

0.5026 

"54 > 0651 

-53,6764 

«53,8i75 

-54.1810 

-53,9992 

0.1818 


o 
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SE-AHT 


LgAST SQUARES DeTeRmINATIOR OF S^IOOTHEP SEA SDREaCB HEISHTS 
CAUCULATED'A’T'IMTERSECTINO POINTS OF ORQSSING REVOLUTIONS" 


calibration AREA 
latitude 59,00 TO 67;00 

longitude 300.00 TO 330,00 

N0‘,0F PARAMETERS ^ i 


NO^ OF SDUTH*T0-N0RTH REVS » $ 

no; of NbRTHSTO-SOJTH REVS *. A 


o\ 



TABLE 13 



NORTh-TO'SOUTH 

Ru>/C‘uJT10NS 

REV NO 

SE3 

6<10 

563 

26'50 

360 

2S29 

479 

4236 

15 

4577 

232 

4562 

?22 

626 

566 

4548 

213 

SOUTH^rO-NOftTH 

revolutions 

REV NO 

'SE3 NO 

640 

563 

2650 

360 

2673 

376 

2829 

479 

4236 

15 

‘ 626 - 

- 556 - 

663 

391 

641 

564 

1167 

101 


o\ 


TABLE 14 



ORIGINAL PAGE IS 
OF POOR QUALITY 


«*»*ih»-ihhhm*CONSTRA IN5D REVS**«*»*«*(Mt* 


R6V NO SeG no 


4236 15 

4i>48 213 


O' 

w 


TYPE OF constraint 
HEISH t SLOPE CURVATURE 

« 

* 


TABLE 15 



A«C data ImpUT FO^ this P'J’J 


REVC/UUTION ^C 


ChOSsiSO 1 

=ori'Ts 

shoothfd sfa 

SURFACE height 

TjHF ^■RDH 

epoch 

nffth. 

-SOUTH 

south- 

•koRth 

LATiTjDE I 

.ONGlTUnE 

nohth-south 

south-north 

RoRth-south 

SOUTH-HORTH 

(FEv/sr(.) 

(KEV/StG) 

tCFGfisSrO 

(HEC.RFES) 

(HtTEnS) 

<KFTERS) 

(HiMUTFSi 

(MIMUTES) 

6-10 

563 

6 4(3 

563 

65,)3433 

346,99063 

22'.4?969 

22,42969 

-1.74C139 

•1.740P9 

16'j3 

36 0 

2656 

360 

05.12945 

325,57243 

6q , Stl56 

60,85156 

8.37898 

6.37898 

2c53 

360 

2829 

479 

65, r|4635 

320,02355 

50.11682 

44,53372 

8.93J82 

14.04771 

<1236 


2629 

479 

65'. n4342 

320.91197 

43,00198 

44,65421 

5,74540 

14,03740 

4236 

15 

4236 

15 

65’, j 2493 

325.36430 

51.92188 

51,92188 

5',2?387 

5.22387 

A’^77 

232 

64o 

563 

04,‘7?®19 

319.71224 

44;B3345 

33,16425 

5'.63083 

-3,22911 

Ab77 

?32 

2629 

479 

64^04519 

323,33159 

52 , 40392 

51,06031 

5.20009 

13.75370 

4562 

222 

640 

563 

03’.».6‘>‘'2 

329,57867 

— 54‘.66899 

56,53514 - 

4,83760 ‘ 

—4,44844 

4562 

222 

2673 

376 

03l 4?44 q 

328,43145 

52'.87i72 

52,48216 

4',9863i 

5'.293il 

4562 

222 

64 

564 

6 1 ■.1334 2 

316,89210 

■ 64'.96653 

41,40091 

6160434 

9 , 03l04 

4 562 

222 

1167 

101 

6li i1®54 

316,83937 

55.67219 

64,53549 

6,61237 

8,643?8 

626 

656 

626 

556 

65^12077 

300,44065 

11^14063 

11,14063 

16',5775o 

16:57750 


■:rv 

(i. 


TABLE 16 



RFVOLUTION NO 
NORTH-SOUTH SOUTH-shfiRTH 
(REv/sesj (REv/sec) 


640 

563 

64o 

563 

2656 

360 

2653 

360 

2658 

360 

2829 

479 

4236 

15 

2829 

479 

4236 

15 

4236 

15 

4577 

232 

640 

563 

4577 

232 

2829 

479 

- 4562 

222 

640 

563 

4562 

222 

2673 

376 

4562 

222 

64i 

564 

4562 

222 

1167 

101 

626 

556 

626 

556 


ARC crossing DATA USED 

CBOSslNG POINTS 
LATiTjOE UONRITUDE 
eOHcReFS) (DBGSEES) 
65 ;j3433 306.9VB83 

- 65’.i2'^45 3?5, 57243 

«5.n'<635 320,82335 

f>I>‘.D'’3'’2 320.91197 

65'.i2493 325,38430 

64,72919 319,71224 

64^94519 323,33159 

- 63^*6482 329,57967" 

63j'4944o 328,43145 

6l, 13342 316,69210 

61^11954 316,83937 

65^i2077 3C0. 44065 


FOR ADJUSTMENTS 

SMOOTHED SEA SURFACE HEIGHT 
NORTH-SOUTH SOUTH-NORTH 
(METERSJ IMETERS) 
22,42969 22,42969 

-- 6o‘.85i96 60185156 

50.11682 44,53872 

43‘.0dl98 44,65421 

51,92188 51,92188 

44.63345 3.3,16425 

52’, 49392 51,06031 

- 54.66899 -- 56,53514 

5?;87i7? 52,48216 

64,98653 4l,4oo9l 

55.67219 64,53549 

il',l4063 11,14063 


TIME from epoch 


north-south 

south-north 

(minutes^ 

(MINUTES) 

-1,74039 

-1.74069 

8.37898 - 

8,37898 

8.93182 

14.04771 

5; 74540 

14.03740 

5',22387 

5,22387 

5;63083 

-3.22911 

5',200S9 

13,75370 

4 ,83760 

-4,44644 

4;9365l 

5,29311 

6160434 

9.03104 

6.61237 

8;64328 

16157750 

16.57750 


cn 



TABLE 17 



y V 


CHARLV-Dy^.AMIC COPf- ALLOCATOR- 

599WORDS OF CoSE HAVE BEEN 
5ISTPj0UTEri among the F0LL0W15}G arbaSs 


SHREV 

8 

KSREV 

6 

sum 

90 

SU^2 

18 

DELTA 

18 

bbias 

10 

kmstr 

2 

ARRAY 

448 


starting addresses oF arrays 

SKREV 127A26000DPO 

NRREV 127436COOO{iO 

SLNI 127A4ACnOOt.iO 

SLH? 127&76COO0riO 

delta 1276200CC0C0 

HFJAS 1276A2COCOC0 

KLSTR 127654000000 

array 127&96C000P0 

arv CHECK aT locattbn 834507 


TABLE 18 



NORTHrTOnSOUTH REVS 


REV NO 

SEG NO 

COEfFIClEMTS 

640 

563 

•0.6583354E Ol 

2650 

360 

0.72303?5E Ol 

4236 

IS 

««0.95955i4E-l» 

" 4377 

232 

0.3065fi33E Ol 

4562 

222 

««0.104495lE 02 

626 

556 

Ol 


SOUTH' 

■to*norTh Revs 

REV No 

SEG NO 

COEFFIClPytTS 

-- 640 

563 

•0,855i 3354E Ol" 

2656 

360 

0.7230326E Ol 

2673 

376 

•0.10639o7E 02 

2829 

479 

0il6522siE Ol 

4236 

13 

0. 

626 

556 

0. 

“ 641 ■ 

■ 564 

»0«34035i3E o2 

1167 

101 

<fO«l5662i6E ol 
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HbVCLUTIONS 
M/S S/M 

crossing PTS, 
Latitude longitude 

(DEGRErS) (DEG3EES) 

PUSgS 

N/9 S/M 

TKETERSJ (MFTERS) 

SEA SURFACE KNIGHTS 
M/S S/N 

(METERS) (METERS) 

SEA SUBFACE 
N/S 

(METEBS) 

HEIGHTS average 
S/N TADJUSTED') 

(HETERS) (METERS) 

residuals 

N/S 

(METERS) 

6<)o 

64fl 

65;i34333 

306,99083? 

-8,5834 

-6,5034 

22,4297 

22,4297 

31.0130 

3ll0130 

3l',0l30 

0, 

2658 

2658 

65,-129446 

325,572426 

7,2303 

7,2303 

60,0516 

60,0516 

53'.6242 

53 I 6212 

53;6212 

0, 

2658 

2029 

651046348 

3201023555 

7,2303 

176522 

50,1168 

44,9387 

42 ,'6065 

42:8865 

42',8565 

Ol 

4236 

2829 

65;o434i6 

32oi?11966 

-O,O0OU 

i;6522 

43,0020 

44,6542 

4310020 

43:0020 

43.0020 

Ol 

4236 

4236 

65 ,'12493i 

325,364300 

-0,0009 

- 

Oi 

51,9219 

51,9219 

51 .'9219 

5i:9219 

51;9219 

0. 

4577 

640 

64:729195 

3l9i7l2238 

3,0053 

-b;5034 

44,6334 

33.1643 

4t;7476 

4l'.7'i76 

41,7476 

0, 

4577 

2829 

64,945139 

323.331593 

3,0053 

1.6522 

52,4939 

5i, 0603 

4974081 

49:4061 

49.4001 

0. 

4562 

640 

63,664015 

329,5786o7 

SIO.4495 

-e';5034 

54,6690 

56,3351 

6571105 

65'. 1185 

65 ', 1185 

Ol 

45«2 

2673 

63,494395 

320', 431454 

(=10,4495 

-10 1 039i 

52,8717 

52,4822 

6373212 

63:3212 

63.3212 

0. 

4562 

641 

61.133418 

316,992101 

sIo.4495 

-34,0351 

64,9865 

41, 4Q09 

75,4360 

75'.436o 

75:4360 

ol 

4562 

1167 

6l. 119545 

316', ’§39371 

910,449V 

»i';'5862 

55,6722 

64,9355 

6671217 

66 : 121 ? 

66:i217 

0, 

626 

626 

65 , 12o769 

300,44:651 

0, 

ol: 

li,l406 

11,1406 

1171406 

4i:i406 

Il’,l4n6 

0, 


0= HeA^ OF RESICUALS 0, 

STA^DARO PcVUTION 0, 


TABLE 20 



REVCLUTIOf<S CROSSING PTS, BIASES SEA SURF I SHTS SEA 3Upi'/(ci'''l!; IQHTS AVERAGE RESIDUALS 


S/N 

N/S 

latitude 

{DFSREES) 

LOf4GITUDE 

(DEGSEES) 

S/N 

(METERS) 

N/S 

(METERS) 

S/N 

(METERS) 

N/S 

(METERS) 

S/N 

(METERS) 

N/S (ADJOsTED) 

(METERS) (METERS) 

S/N 

(METERS) 

6A 0 

640 

65,134333 

3Q6,?98833 

-8)5834 

-8', 5834 

22.4297 

22,4297 

31,Cl3o 

3l’.0130 

3lloi3o 

0 , 

6Ao 

4577 

64,729195 

319V712239 

-8,5834 

3,0658 

33,1643 

44,8334 

41.7476 

41^7476 

41’,7476 

0 , 

64 g 

4562 

63,664815 

329,578667 

-8,5834 

-10:4493 

56,5351 

34,6690 

65, '1185 

65', 1185 

65,1185 

0 , 

2658 

2658 

65,129448 

325?572426 

7,2303 

7:2303 

6tf,S5l6 

60i85l6 

53V6212 

33'.6212 

53’.6212 


2673 

4562 

6s;494395 

5?8;«31454 

b10,839i 

“10, "4495 

- 52.4822 

32,8717 

63.3212 

6313212 

6S;3212 

0 , 

2629 

2658 

65;Q46348 

320SS23555 

ii 65 Z 2 

7:*23o3 

44,5387 

So, 1166 

42.6865 

42'.8865 

42.8865 

0 , 

2829 

4236 

__ 65,043416 

320. 911968 

i' 6523 ^ 


44,6542 

_ 43,0020 

43 ‘,’ 0 i)ao 

43', 0020 

43‘,00?0 

...O'. 

2629 

4577 

64,945169 

323.331593 

1,6522 

3:o85« 

5l, 0*03 

52,4939 

4974081 

49’.408i 

49,4081 

0 ) 

4236 

4236 

65;i24931 

325,384300 

0, 

-o:oooo 

51,9219 

91,9219 

Si. -9219 

Si*.9219 

5i'.9219 

0, 

626 

, 626 

65.120769 

3009440651 

Oi 

o: 

11,1406 

lt,:406_ 

li',‘1406 

ll'.14((6 

11’. 1406 

0, 

64i 

4562 

61.133413 

316,892101 

s34 , 0551 

910^4499 

41,4009 

64,9865 

75 , ’4360 

7514360 

75‘. 4360 

0, 

1167 

4562 

6i.'110545 

316,839371 

-1,5862 

-1074493 

64,5335 

95;6722 

60;i2i7 

66*,1217 

6&;i217 

0, 


<£> 



TABLE 21 



FORMAT OF SEAHT OUTPUT TAPE 


[EADER record 


WORD NO . FORMAT 


1 

2 

3 

4 

5 


I 

I 

I 

I 

I 


DESCRIPTION 

Dummy ( 0 ) 

Dummy ( 0 ) 

No. o£ Data Records 
Segment No. 
Revolution No. 


>ATA RECORD 


WORD NO . FORMAT 


1 

2 

3 

4 

5 


I 

I 

I 

I 

I 


DESCRIPTION 


Time Ci’^illiseconds) 

Latitude (milliseconds) 

Longitude (milliseconds) 

Sea Surface Height (micrometers) 
Dummy (0) 


TABLE 22 
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3.6 


LIMITATIONS 


Several dimensions have been set in SEAHT to values 
sufficiently high to accommodate all runs that are expected 
to be made by the program. The dimensions set produce only 
the following limitations: 

• Number of North- South passes < 200 

• Number of South-North passes ^ 200 

• Number of data records per pass £ 2000 
(—33 minutes at the 1/second data rate) 

• Number of Selects/Deletes from one data tape 
< 50 

• Total number of input data tapes ^ 6, 

Within limits, all the above dimensions could be significantly 
increased if needed, without increasing the total core used 
by the program. The arrays used in the parameter estimation 
portion of the program are all dynamically allocated, and so 
the number of passes that can be processed depends upon the 
number of parameters estimated per pass, and the amount of 
core that is available. 

SEAHT will not function properly if: 

• Some revolutions do not cross either a constrained 
revolution, or a revolution which crosses a 
constrained revolution. 


71 



rhere are no error terminations in SEAHT, and improperly setup 
jobs may be terminated either via system abort, or by meaningless 
program output. 
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SECTION 4.0 
PROGRAM LISTING 


A listing o£ the working version o£ the program follows. 
The placement of LINK cards is shown. 


73 
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1 

2 

3 

4 

5 

6 
7 
a 


11 

12 
13 
11 

15 

16 

17 

18 
1 ? 
20 
21 
22 

23 

24 

25 

26 
27 
2§ 
2? 

30 

31 

32 

33 

34 

35 

36 
3/ 
38 


CSEAHT- 


DIMENSJOH IBOREiSj. 

COMMON/TAPEg/lN|TP;i0UT,C3TAP, MSlAPEf SfjTAPE;sCRB*SCRC,DATAPE,DOUT 

COMMON/KONTBU/NOPaRi NOM3,NOSN, IKf^STR,^EOIT, max 

EQU I valence HCOHEd) , ICSN) *( I qORi^ ( 2 ) * I CNS ) , ( IC0R6<3) ,NSUM1) , 

. ( ICQRfa{4 j J VSUMg) / ( rcOR£(5 > , NdeU > flCORe ^6 ) , NP) , 

, ( IC0RE(7) ,NKNS TR J ( ICD ^E^ 8 ) , I MAX ) 

DATA nQcRE/6/i ImOrE/0/ 

CALL LLINK(5HLINK1 J 
CALL SEARCH' 


A!P=NOPAR«(NaSiV + i) 

MSUMl=NP«(Nfi+li 

^|SUM2Bf^p^>2 

NDEL?NSJM2 

ICSNtNOSN 

ICNStNONS 

NKNSTR=lKNSTRtt2 

IMAX?MAX«4 

DO IS lal.NOdRE 

IHORE=IMORE!JiICORE( p 

10 IF(MODUCORB{ I P2) ,NE,0> J CqrE U ) t J CRRE ( I ) +1 
KOREB 3320 

write ( 6, 9) IMORE/ (IC0RE(K)iKall8) 

5 FORMaT< flpajxPcHARLYcDYNANlC CORE ALLOCATOR-^ P //, 16 , 

pWORDS OF cgSE HAVE BEEN * , / P D I S TR IBUJED AMONG THE FOLLOWINQf, 

. , ’\>ARRAX5p/^;iOX, »SNR=V ' 1 5X / 1 61/ ,ilOX » NSREV » » 5X 1 1 6 , / , lOX , JSUMl' , 

.6X*l6P,iOXpS'JM2i , 6X; 161/ilOX#^ 'PELTAP5X, l6iVi lOXPSBI AS ' ,5X, l6. 
./,10XPXNSTS> PX, 16,/, lOXi 'ARRAY S5X„I6) 

CALL LLlNK(5HLtNK2 } 

CALL LLINK(5HLINK3 ) 

CALL CdARLYfKOREB, I MORE, NCORE , I CqRE P'OREG , $100 ) 

GO TO log 
100 WRITE(5,101} 

ni formaT( fp , tqHaRLy has been uNsyccEsgrUL in alloca'^'Ing core,') 

102 call LLINK(5hLINK4 ) 

CALL solve 

STOP 

END 
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1 

2 

3 

4 

5 

6 

7 

6 

? 

ly 

11 


CBDATA 


block Data 
INT eQER (?TAB^DnUT 

COMMON/KONTRL/NPARCiNDNSiNOSN, tKNSTR^60|^;,MAX 
COMMQN/TAP68 /INTP;iOUT’,GTaP, ISQRt5)#R0UT ,, 

DATA InTP/ 5^;IC)UT/6/» l5CR/Uil2,a}3,lgi,,i5i',GTAp/10/*DOUT/lO/ 
data NBARC/g/iMDNS/0/,MOSN/Q/ 

DATA IKNSTR-ifO/ 
data E01T/1?IQ,/ 
data MAX^O/ 

END 
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1 

2 

3 

4 

5 

6 
7 

a 

? 

18 

11 

12 


CNUMBER 

INTEQER RUNKTION NUM05R(K, lAi ID) 

dimension latiDJ 

NUMBER?0 

IFUD.LTar RETURN 
DO 10 

1F<K,EQ, U(’>n GO TO 20 
to CONTINUE 
return 

20 NUMBER*! 

return 

END 


a\ 



s ^ 



ORK CTNKT 

f 


\ 


OOQODDQIlDOOOIlODillillOOtlOOOIlDOIlOIlOQDOOOOIlODflOSODO&ODOHMIlilOO&DODailOIIGaOOOOOOODDOOO 
I ] 1 t I t I I t It II u I) II IS iMMi II II II I) II V IS ’s II M II II 31 II a H I! i( II 11 a ii ii it ii m <5 II II II II SI si si m h ss « ii a si u ii ii » ii is ii ii ii ii it ii ii ii ii is n ii ii ti it 

11 1 1 111 1 1 1 1 1 1 1 1 M n"n 1 n 11 ! 1 1 n n 1 M 1 n t n ! 11 m 1 n 1 1 1 1 n n 1 11 1 11 n 1 11 n 1 11 1 1 M 

2222222222“22n2n"22J222222222222222222222222222222222222222Z2222222222222?2222 

"323333''3333333~32333323333333]3333323a3333]33333333333333333233333333233a333333 
4444444444444444444144444U444444444444444444H44U44444444444444444444444444444 
S$S9555S9~55S55S5~5555555555555S5S5SSS5SS5555SS55S5355553553S535S55S55SSSSS5S5SS 
6SS6G6666GG86G3BCG6B&GSGEGGS6GG06G66E66CCGGSSBGG6EG6EB8GE5SS(G66GES8E6{EeE6GE663 
7 7 2 7 7 n n n 7 7 7 7 7 7 7 7 7 7 n 7 7 7 7 1 n m 7 n 7 7 7 7 7 1 n T 7 n n 7 7 7 7 7 1 7 n n n 7 7 7 1 1 7 n m 7 n 7 7 7 


9 9 B 9 B S S 9 ~ 9 9 9 9 9 9 E ~ 9 9 9 9 9 3 9 9 9 3 9 3 9 3 9 9 9 9 9 9 9 3 9 S 9 3 3 9 9 9 9 B 8 9 9 S 9 9 9 3 3 3 9 a 9 9 B 9 B B 9 3 9 9 9 9 3 3 9 9 9 9 9 
I 1 3 4 s t > 1 I It II II II ms IS It II H It II II II II IS » II IS » II II II 11 n IS II II II ’I It II I! II 11 IS It I' II IS sc SI SI 'I SI IS SI SI SI SI It i: ii ii ii ss ii si ii it ii n ii ii ii is it ii ii n ii 
NEcr-ioai — ■ - . - .. ^ — ,, , , — 


57355 01 U 3 - 16-77 16,999 


1 

t 

6 

5 

5 

6 
7 

9 




1 

U 

12 

J.3 

14 

13 

16 

1/ 


IS 

1 ? 

20 

21 

22 

23 

21 

25 

26 
27 


29 
2 ? 

30 

31 

32 

33 

34 
3? 

36 

37 

3 5 
3? 

40 

41 

42 


qsearch 

subroutine search 
external USBRR ■ 

LOGICAL MINUTE 
LOGICAL NOPATA 

integer 0ATi((8J ,SnTaPE,ScR8,REVNQ,CDNaMEU 1) i'CAR0lD,TAPLlM,UNIO 
integer SNRBV<200 ) iNSREV< 200) ,CRSNGSf5Q0Q,2h'MONLEQ(200 »2)i 


t 


I 


NDELl(2)i:NDELT(2) I lNTNOS{20QI2J 

integer patare 

integer SELgbT(50;6 ? , delete ( 50*6) iKN0t6 5 IkNA{6J , ITaPE(6,2) 

REAL MnUAT7.MAXLAT,MnLON,MAXLQN 
dimension WaRDS(4),DTl(4 ) ,SShT(2) 

DIMENSION K5STR(5000i 2) 

COMMON/'kOnT8l/NOPAR#N 3NS,NQSN, iKNSTRe-GDIT.HAX 
COMMON/ALT/Bl,iRRAV(4,lOOO) ,Nj;ARRY<g,.10 00 ) 

COMHON/TAPE$/INTP;ioUT,GTAB,NSTAPE,SNTAPE,SCRB,SCRC,DATAPE,DOUT 
equivalence tXcR0SS<DATA<3)) , (YcR0SStDATA(4n , (QJ,DATA(6) ) » 
(TIfOATA{7>), (TJ|DATA(8M, (6UDATAt5)i;sSHT(l)) 
equivalence (UnIQiW0RDS(4) > 

Data NQDATA>i, false,/ 

DATA MINUTE^iTRUe,/ 
data MINUS/;!/ 

DATA NDELT/§»»0/.NCARD/ii/,TAPHH/6/,T;TPl/0/» JTAPE/12«0/ 

data CDMAHE/6HaREa ,5HPaRAM , 6HN0DAiT,A i 6HcONSTR , 6HED I T >6HDATA 

6HSELEeT;6HDEU&TE;6HVERTCLj.6H2ER§ ,6HMINUT£/ 

data KNO/6»?/,KNA/6ttO/ 

data MINLAT/-90 .0/ 


MAXLAT/ 90 ',/. 
MINLON/COO./; 
MAXLON/ 360 ,/ 
data Z6rq/o;q/ 

DATA IKMSTR/d/ 

DATA EDITl/ 1000 ./ 
INTEGER R,n 


20 RFAD(,lNTR,lfQ 2 ,PND» 2 l) CARD I D , IARRAY/.DT 1 
1002 F 0 RMAT(A 6 , I!J; 4 E 15 , 5 ) 

I 3 NUMBER(CABDID,CDNAM 5 ,NCARD) 

IFU .NE, 6 > Go TO 20 

itpi*itoi+i 
iTAPEUTPiaielARRAY 
K = DU( 1 )+ 0,1 
ITAPE(ITRli 2 )=MAX 0 (i<,l) 


I 
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'••4 

O 


43 


46 

47 

49 
4? 

50 

51 

52 

53 

54 

55 

56 

5 7 
58 
5? 
60 
61 
62 
6^ 
64 

63 

66 

67 


69 

6 ? 

7Q 

71 

72 

73 


76 

77 

79 

79 

80 
as 
aa 

83 

Q4 


21 

6 


1 


2 

3 


4 


5 

7 

B 

9 

10 

1003 

11 

50 


GO TO 20 
REWIND INTP 

READtINTR4^02VeNDs50) CARDID. U8RAY«0T1 
I«NUMBSR{CA8DIPiCDNAM£’,NCARD) 

GO TO Ui2i3;4',5i6',7,8,9,ig,il)/I" 
MINLAT?DT1(1) 

HAXl,ATsDTl(2} 

f4lNLON:;DTl{3) 

MAXUQN=nTl(S) 

00 TO 6 
N0PAR=IARRA? 

GO TO 6 

nodata?,true; 

ifiiarraX.gt^o) scrb^iarray 

READISPRB) 8OMS,N0GN 


Go TO 6 

iknstrsiknstr+j 

KNSTRUKNSTB;i)sDTl(l)<^iq0y00,O + DTl<?)«l0,0 
KNSTR(IKMSTfi;2>=lARRAT 
GO TO 6 
EDIT=DT1U) 

GO To 6 

call SELDEUhTPl.lTAPc!, IARRaY^&TI^SEL^BcTXKNO) 

GO TO 6 

call S6LRELHTP1, ITAPS, URRaY, PTl, DEUBTE^KNa) 

GO TO 6 
ISETVpUI 
GO TO 6 
ZEROtDTUI) 

WRIT6(6,100?) ZERO 
F0RMATUX,»A zero option used WiTh'iR 13,5, *L,ONqlTUDE 2ERO-GET»> 
GO TO 6 

MINUTE*, FALSE. 

GO TO 6 

IF(NODATA) RSTljRN 
00 55 Ul.ITRl 
lUNlTsITAPElra) 

NTPASSalTAPBU.2) 

CALL FLGPRCtlUMlT/uSERR) 

DO 90 Ksl.NTPASS 

reao(iunit,snd>»90) words, revno 

REVNQsREVNO^JlOnoOO + UNIQttlO 
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H6 

87 

89 
S? 

90 

91 

n 

93 

94 

95 

96 

97 
99 
9? 

i0i3 

101 

102 

1'03 

1U4 

lo? 

10^ 

IQ? 

10? 

10? 

ue 

lu 

U2 

u? 

u? 

11? 

116 

11? 

118 

11? 

120 

12| 

12? 

X2i 

12f 

12 ? 

126 


16.99? 


IF(NUMa 5 R(RBVNO-SELECT(liI),KNQU>),fiT, 0 )f GO TO 60 
IF(KN0{U ,Na‘,05 GO TO 61 

IF(NUMBEfl(R0VNO,DELET5<l,n,KNArt) ) .iO.OV GO TO 60 
61 READ( IUmTiGNDa90 ) WORDS 
QO TO 61 
60 NRsQ 
iPRTUsO 


H=1 

65 ReAUUUNlT,0ND»78) WORDS 

IF(W0RDSl2)^LT.^1INLAT.^R,W0RDS(2) .GTO-HaXLATJ QO 
lF{WORDSr3)*LT.MINLON,OR,WOROSl3 ) .qT^MAXLON) QU 
IF(MINUT0? WORDS {1)=W0RDSU)«144Q, 

ARRAY 

ARRAY(2jH)=W0R)3S(2> 

ARRAY{3iM)?WdRDS(3> 

ARRAY(4,M)pW0R11S{4) 

IF(M,L,G,2) eo TO 74 

IFU)IFFh>UR$AY<2iM)«ARRAY(2,'1-1)),GE90,0} GO TO 73 
NR*NR+1 

BACKSPACE lONir 
IPRTL=1 
GO TO 79 

74 IF<M,EQ,3) DlFrsARRAY<2,2)^ARRAY(2,l| 

75 NRsNR-H 
MsM + l 


TO 65 
TO 65 


IF(M^U6,100$) GO TO 65 
73 IF(NR,l,E..0) Oo TO 90 
79 IF(NR.Ue.l.AND, 1PRTI,3T.0) GO TO 60 
IF(NR.GT,HA«J MAXbNR 

1F< I SGTVreOSlU, And, (ARRAY (SiGJ^'aRRAlli 3,1) ) ,QT,0i 0000001) GO 'TO 80 
X5sAfi!RAY(2,?|)-ARRAY(2ri) 

IF(XS,GT*0,g) go TO 80 
NM=NM^4 

NSREy(NM)3RgyNO 
WRIT5<N5 TaPS) ReVNOiNR, array 
IF(IPRTu„QT 5(3) GO TO 6g 
QO TO 90 
ao NN»NN+1 

SNREV(NN|sRBVNO 
WRITe(S'MTAPB) R£VnO»N4, array 
IF( If?RTL,GT5d) GO TO 6o 
90 CONTINUE 
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%27 

m 

X2^ 

130 

131 
13? 
133 
43i 
133 
13<i 
S37 
13a 
13? 

140 

141 

142 

143 

144 
143 

146 

147 
140 

14? 

150 

151 
15H 
153 
«5^ 
153 
15g 
157 
139 
15? 
160 
161 
168 

163 

164 
163 
166 
167 
169 


55 continue 
eNDFlLg NSTAPE 
ENDFILB SNThPE 
REWIND N3TAee 

WRITE(Sc«B) NM.'nNi’NSRcVjSNREV 
DO 280 

REAU(NSTApE?,ENOa200 ) REVNOfN I ; ARpAY 
lF<2aRQ,BQ,9;cn go to 56 ' 

DO 57 Ual,N*I 

if((array<37ui-zfrd> , 3a;o.g) qq to sg 

ARRAn3,U)=ARRAY(3.L)+360,0-2ERO 
GO TO 57 

53 ARRAY(3,U)aARRAY(3it.W2ER0 
57 CONTINUE 

56 rewind SNTAGfB 

DATA(2)«N 
EXUNGMsaRRAJO.'X) 

EXUATM? ARRAttgjlJ 
ENLNGM?ARRAY|3 ,nI ) 

ENlATMTARRATt2.NI) 

DO 150 N*1,SN 

READ(SNTAPE;ENI) = 150) REVNQiNsI.ARRY 
IF(ZERQ.EQ,0'.O) Go TO 59 

DO 70 L«1#NU ^ 

IF<<ARRYt3,fc)-7ER0) ,G=,0iO> GO TQ 69 

ARRY(3i*l.)aABSYI3»L.>*36 0,0-*ZERO 
GO TO 70 

69 ARRYt3f-L»BABRY(3iL)-ZSR0 

70 CONTINUE ^ 

59 DATAUMN 

ENl,AIN?A8RYI2il) 

ENLNGNrARRYlSil) 

EXLATN?ARRY|2,.MJ) 

EXUNaN*ARRYigiNJ) 

D7=(£NLN«M-gNLNGN)o(EXUNGMHEXUNGNj 
D8a(EXLATN-eMLATMH{EUATN''EXl,ATNl 
IF(D7*GT*.0,9,OR.D8,GTiO,g) GO TO 150 
CALI, CROSS (>[CR0SS,'YCR05S7ENLATMiENUN(^M,EXLATM.EXl.NqM,ENLATN, 
ENLNgN|ERl.ATN,EXLN3N, INDC) 

IFUNDC.EQ.I) go to 150 

CAUL GHTlM(8CRDSS;yCRD5S,Gl,GJ,Tl.Tj) 

1F<Z£RQ,EQ,(3:0) GO TO 72 
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116 ? 

173 

171 

172 

173 

174 
173 
17^ 
177 
179 
17? 
109 
161 
182 

183 

184 
103 
186 


18 ? 

199 

191 

xn 

193 

*94 

193 

196 

197 

198 
19? 
? 0 « 
201 

m 

?03 

?04 
303 
306 
30 7 

3oa 

?o? 

2ifl 


IF< acfio3S+?ERO),GT, 360,0) GO TO 73 

YCROSS¥yOROSS+Z6Ro 

GO TQ 72 

73 YCROSSaYOROSS~360;0+2=RO 

72 IF(ABS(S 8 HTii)),GT.EDITi) GO TO 150 
RGSlDsARS(S3l^T(2>»SSHrfi) ) 

IF(R 6 Sln,GT'JED!T) GO TO IgO 
WRITE(SCRB) data 
NTATA sNDATA^l 
CRSNGS(MPATA;1J3 DaTA<2) 

CRSNGS(NlDAT'A;2)=DATAa) 

150 CONTINUE 
DATaU) ant 
VIRITE(SCRB) data 
200 CONTINUE 

ENDFILE 9CR8 
WRlTE(6,200?J 

2000 FORMAT (»lt, ei3)(, »SEAHT >,////, 36X^ 

i>LEAST SQUA 8 ES DETERMINATION OF SMOOTHED SEA SURFACE HEIGHTS’,//,. 
237X, ’CAlGUlATEO AT INTERSECTING RoINTS OF; CROSSUG REVqLOT { 0 NS» ) 
WRITE ( 6 , 300?) MlNLAT,MaLAT,MlNLgN»MAXLON 

3000 FORMAT(////y57X, •CALI3RATION AREAl » > LAT I TUEF. f,F8,2| 

1 » T0»,tF8,2j:/,5lx;tL0NGITUDE UF8,2#1 TQ»»F8,2) 

WRlTE(6,300f ) NOPAR 

3091 FQRMAT< »0» ,?3X, ’NO.OF “PARAMETERS a i/H,) 

WRITE ( 6 , 3003) NN 

3003 FORMAT(.///,50X, »N 0 , O’ SOUTH'*TO-NORTI:I REVS s 1,12) 

V|RIT 8 ( 6,300?) NM 

3002 FORMAT! »Qi ,^9X, »NO, O’ NORTH-TQ-SOUTI;) REVS « »#12) 

WRITE! lOUTitOO'j) 

DO 67 lalfNM 
INS=NSR5V! U/IOOOOO 
lNSS6Gi={N3RgV! p»iNSttltl000 O)/10 

67 WRITE! IOUT,iOO$) INSiINSSEG 
WRITE! IOUT#$007) 

DO 68 Isl,NCJ 

ISNsSNREV! n/lOOOOO 

ISNSEG = {SNREV( I )iilSN^aPOOOO)/lP 

68 WRITE! nUTil'OOi) iSNjISNSRG 

1005 FORMAT! ’!» i50 X, ’NORTH«To»SOUTH REVOLUT IONS ’ t / , 54X , 

, ’REV NO SEG NOP /) 

1006 FORMAT!’ r , 54 X , I S5 , 6 X , 1 5 ) 
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00 

oi 


211 


3if 

?l3 

314 


w 

31^ 

317 

319 

31 ? 



323 

323 

324 
329 
326 



32? 

338 

331 

338 

333 


334 

3355 


336 

337 
339 
33? 
343 
341 
348 
343 



346 

347 
340 
34? 

350 

351 


353 


1007 F0RMAT( «Q> i50X.'’S0UTH«TOsNORTH RaVQl,UTJONS» ,/,^'4X| 

,»R6V NO . SSG 

V^RlT6<6,399?} 

DO 3007 K3lj;jKNSTR 
KNREV = KN3TRIK^U/100000 
KNSEQs<KNSTPi<i<a)*KNRSV«10QOOO)/lO 
IF(KNSTR(K*2> ,RQ,1> W4 6, 30Q4.) KNRBV,KNSEQ 
IF(KNSTRIKi 2> »EQ,2> WR I JE < 6, 3005} KN6EVjKNSEG 
IF(KNSTR(K,3) ,ga,3) WRITE{6,30P6} ,KNgEV»KNSEG " 

3007 CONTINUE 

2999 F0RMAT(*1'i 67X» if CONSTRAINED REVSft*'tHHMMHMHH»**//’/>43X# 
,«REV N0»;5XJ:»SEG NOMOX,«TYPE of COySTRAlNTi, 7, 65X, ‘HEIGHT* ,3X, 
,tSL0PE*,3X,rCURVATUREI,/7> 

3006 FORMAT (43X,iI 3, $X, 1 5» lOX, U ' , 7Xi * 9^ , 98, » ) 

3005 F0RMAT<43X,55,AX, I5,10X, itt* ,7X;*4r) 

3004 format ( 43X,T5,6x»l5|10X, *#» ) 

WRITE(IOUT,7000> 

210 WRITEUOUT»7001) 

WRITEUOUT,I002) 

WRlTEUOUT.roOH) 

REWIND 3QR0 
READCSCR8) -n 

220 READ(ScRa.EGld = 240) DaTa 

IF(DATAU) ,E0, MINUS) GO TO 220 
U = DATA{3) 

12=DATA(1> 

DATAU)sNSR0V(I1) 

DATAl2)sSNRBVn2i 
INSaDATAll) 710001)0 

INSSEG?(PATAU)-INS<»100000 )/10 
ISNbPATA<2)7100o 00 
ISNSBG^{DATft^2v-ISNftlOOO{)0)/10 

WRITEUOUTiI0 0 3 )INS» INSSEGi.ISn;I9NSE[^..DATa(3),0aTaM)»DATA(5), 
,DATA<6),0ATJ!(<7),DaTA(8) 

GO TQ 32P 

240 NPRiNTtNRRiSt+i 

IF(NRRINT.GTVI) RETURN 
300 ndeludso 
NDEUI (2)?0 
11 = 1 
1122 
NUIMtNM 


o S 
nJ S 

t?3 




CO 
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•(» 


?53 


?53 

?5<? 

357 

259 

35? 

360 

361 
362 
363 

w 


365 

366 

367 
263 


36? 
379 
“ 37^ 
373 
373 


374 

379 

376 

377 
373 
37? 
389 


381 


38 8 
389 
33-< 
389 
386 
38? 
398 
38? 

ii. tt ^ H 


399 

391 

393 


330 DO 650 HliC^LH 

IF(NUM 8 !=Rn;:NOvJLEQ(l,iniNDElTU'U),WE ,01 GO TO 350 
1457 DO LOOP INDEX I HaY NOT BE IEDEPINED CAuL OR ABNORf'AU- TUNCT JON 
C DETERMINE THE NDMBER OF C= 50 SSINGS QF' ITH NS(Sn) aRO 

ncross^o 

DO 340 J*l, 8 flATA 
IF(CRSNG 3 Un n . NE, n 00 TO 340 
C TEST TO SEE IF gN(vlS) ,ARC DELETED. J . 

JF(NUMBER(qRSN>' 5 S( J, ID ,N 0 NLEG< 1 ,U>jNDELT(IU) ,EC. 0 ) 
,NCR 0 SSsN 0 R 0 SS +1 

340 IF(NeRQSS,GE‘,NTPAR 5 GO TO 350 
IFINCRQSS,E 0 .QJ go to 345 
1 T=NSREV( 1 ) 

irUT, 6 Q. 3 > ITiSNREVU) 

ITsNWMB£R(IT;K’JSTR, IKNSTR) 

IF(IT.EO,> 0 ) go to 345 
IT = INCR 0 SS*D-'<NSTR{ |T, 2 ) 

IFdT.GE.O) GO TO 350 
345 NDELI(in=NDELtnn*l 
NDEETUUsNDELTC ln +1 
J = NDEUTan 
NONLEGtJDlUl 
350 CONTINUE 
IlslUl 
11 = 1 
NLIM=NN 

IF(II,LE, 2 ) flO TO 330 

lF(NDELIU)J!flT,D, 0 R,NDELI( 2 ),GT,!J) GQ TO 300 
l 2 aMAX 0 {NH,HN) 

DO 410 I? 1 ,T 2 
DO 410 J=il ,3 
410 iNTNQSlIiJ) s I 
NLIMtNM 
DO 420 Hsl ^2 
ILlMrNDEUmi ) 

IFUUM.BQ.fl) DO TO 420 
DO 415 I? 1 ,SLIM 

415 IF(NUMBERU;NiONLEGa,n).Il,IM),Nt^, 0 ) INTNOSUiIDbO 
1457 DO LOOP INDEX I MAY nQT BE REDEFINED iN CALL OR ABNORMAL FUNCTION 
420 NLIHtNN 

DO 430 Ifll.NM 

IFUNTNOSUn) .&Q.Q) ^0 TO 430 
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oo 

U1 


393 

?.9% 

396 

397 
39^ 
39? 

|oa 

9oi 

303 
?Cl3 

304 
309 

306 

307 
30? 
309 
31B 
311 


31? 

313 

314 
3l9 


316 

317 

318 

319 
330 


331 

323 


525 

324 

329 


323 

32? 

330 


331 

332 

333 


334 


nonstnoms+i 

NSREV{N0NS)^NSr^e9(l) 

INTNQS(NONSa>nINTNOS('l/i) 

430 C0NT1NU6 

DO 435 I»1,HN 

IF(tNTNOS<l72).eQ,0) SO TO 435 
N0SN?NQSN+1 

SNREy(NOSN)iiSNReV{n .. 
IMTMQStNQSN#:2>«INTN0SU »2) 

435 PONTINUS 
NOKNST=0 

DO 490 |*=1,1KNHTR 

KaKMSTRUai 

IJ=NyM0SR<KXNSReV;NOMSJ 

IF(N,EQ,'0) lJ«NUM8ER<K,SNREV,NbSN) 

1F(IU,EQ,.0) GO TO 480 

NOKNSTsNjOKNST + l 

DO 475 J*1,S 

475 KNSTR(NOKN'S.T;j>aKNSTR( a J) 

480 CONTINUE 


IKNSTRfcNjpKNST 
REWIND SORB 

WRlTE(DATi^Pe) moNSVnOSNm (NSReVU) , I =1«N0NS ) , ( SnReV { I ) i I=l»NOS|y}> 
, MKNSTRU.Uli iPliNOKNSnjjnia) 


READ1SCR8) 

510 MDATAag 

550 READ(SCR0,ERDs6OO) DATA 

IF<DATa( 1> .go, MINUS) SC TO 590 
DATA(l)aNUM0ER(DATAU), INTNOS(i;2),N0SN) 
DATA<2)aNUMBER{DATA(2) , lNTNOS,NONS) 
IF(A8S<S3HTU)) .GT.EDrU) Go TQ 350 
NDATAsNOATA?! 

write(datape) data 

GO TQ 150 

590 IF(NDATArGTlO) WRITEC DAT APE) DATA 
00 TO 510 

600 ENDFIL8 0ATAPE 
SCRBtDATAPE 



t 


WRITE(ICUT# 700*3 ) 

GO TO 210 

6000 FORMAT! fl»,l5 X,' 'The CONSTRAINED ARCUl7aHAS BEEN DELETED’) 
7000 FORMAT(iHl,3DX, »ARC CRCSSING DATA INgUT FOR THIS RUN’) 
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533 7001 FORMAK »gi ,7x, *REV0L,UT10N NO M2X , I CgOSS I NQ PO!NTSt,6Xi 

536 1‘SMOOTHHC sea surface HEtGHT»,?X«lTltlS FROM EP0CH»>, 

537 70D2 F0RmAT(» ♦ »8X, tN0RTH»30UTH " SOUTynNOgTH* ,5X, f LATITUDE LONGITUDE’, 

538 24X,.»N0RTH-SlJiJTH SOUTiJ^NORTh’ NO gTH-SoUTH SCUTH-NORT'H ’ > 

33? 7003 F0RMAT<» »,l5,2Xn5,’ M 5«3X , FlO , 5 , IX » PI 0 ,5 , 3X , FIQ , 5 , 

540 ,3X«.Fl0.3j;5X^F10.5;4XiriO,5) 

341 7004 FORMAT(tHl,30X, ' ARC CROSSING DAT, A USED FOR ADJUSTMENTS’) 

54? 7005 FoRMAT(J * i 3X» REV/B = G )>*, 4X, » <R^V/SgG >' T6X, t (DEGREES) MX, 

343 , M degrees ) f'fiX, '< METERS ) ^, 5X ;•< METER? P , 8X , M MINUTES )» ,4 X-, 

344 .MMINUTES)’} 

349 end 

1470 equality QR NON-EDUaLITY comparison May n«t be meaningful in logical if 
7 MEMORY EXPANDED, USE $LIMITS QR CORSa OPTION FOR NEXT RUN 


CO 

c^ 


EXPRESSIONS 
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03 


i 

‘i 

i 

4 

i 

6 

7 

§ 

? 

la 

It 

4 ? 

14 

l 4 

19 

1 

1 

19 

1 ? 

2 a 

2 1 

22 


SU8RQUTIN6 SeL02t<N0TPSj!TAPS,NaRRAX^S£tEcT,MS€i.} 

INTEGER £|ELB£;T{ 90 ;i>'VlTAPS<i),NS 6 l,U^ 

OIWENSION A^RAYtl) 

IST=4 

lUMfNOTRS' 

iF(N.ea,g) to 3o 

lST«NUMRgR<S^lt’APS»MOrPS» 

IF<IST,SQ,0) RETURN 
35 ILIM’IST 
30 DO 50 IT?IST,;iLI« 

<»NseuciT) 

DO 40 

lS6UTARRAY(y)»loOOOO,0*ARRAY(J>tl>«lO^C! + Oti 

iFUSeU^EQ,?) SO TO 40 

K»K *1 

56USCT<i<ini4lSEU 
40 20NT3NUS 
NS6U<ITHK 
50 CQNTINUe 
return 

END 


§1 

top 

2 

e tKl 

as 



Oi 

z 

3 

7 

§ 

U 

. i2 
13 

15 

16 
17 
IB 

oo 1 ? 

OO 20 
21 
22 
23 

25 

26 
27 
26 
2? 

39 

31 

32 

33 

34 

35 

36 

37 

38 
3? 

41 

42 


03-16-77 

CCSOSS 


17.009 


.L-NUNQlN)yE?iLA.TM,E)^LNGM, INnc) 

nOHMON7Al,T/ MI/ARRAV {4,iooO>iN>J» ARRYJS4,1Q00> 

HI^LATMfAMLATM 

f!MLNCH*A^{|LNSM 

I:XLATH?A!VCATM 

r.XUNQM5ANLN(5M 

IMDCTO 

IF|BXLATM,ER;ENl.ATN,AND.EXL,NGH,EQ.ENl,hlGN> GO TO 501 
ir(EXlATN,e®.eNl.ATM,AM0»eXLN3NtE9,ENl;NGft) QO TO 500 
^FtEMUATM.LTiexi ATW? SO TO 10 
ao TO 20 

5 IF(fcXLArM,G,r;'EHl ATN) 30 TO 30 
C,Q TO 40 

10 liQ 11 Ial,N0-l 

11 JF{ENLATM.Gg;ARFV(2« I) ,AND»EMLAT«.Ue'i’ARRY(2, I+i; ? 

•OO TO 15 

GO TO 99 

I£:^^°^^W‘-A^M•&^’R«<'<'H^ARRY(2,^^^)lARRY(3,^.ARRY{3^J*l}.Y^ 
ir<Y.GT,6Nt.SGM) GO TO 99 
CO TO S 

20 t'O 21 lal^NT.l 

21 )F<EXLATN.l,0,ARPAy^2, 1) ,AN0,EXUATN,GB^ARRAY(2 i I+lM 

.00 TO 25 
GO TO 99 

25 CALL YFrND(EXLATN;ARRAV<2, i) URRAYU*! *1) , ARRAY(3, I ) > ARRAY(3fI*l) , 

IF(Y.UT,EXLSQN) GO TO 99 
GO TO 9 

30 DO 31 I»l,Ny«l 

31 IF‘eXUTN,Qg;ARPY(2jn,AiND,EXLATf1,UE?ARRY(2|I + m QO TO 35 
GO TO 99 

)'E;l5'5^S^l-'^^^'‘A'^'^''<2tI»»ARRYUfm)iARRY{3jniARRY(3iI+U,Y) 
IF<YtLT,£iXL8QM) GO TO 99 
GO TO 6 

40 DO 41 Jiil.N'Iel 

41 IF(ENLATN.|,E.aRRAy<2i U,AN0,gNLATN,GE';ARRAY(2, I+IJ) GO TO 45 
GO TO 99 

45 CALL- YnND(E«LATN;ARRAV<2,n,ARKAY(2*iai.ARRAY(3,n,ARRAY(3n*l)» 

' IFCY.GT.ENUlitaN) Go TO 99 
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057 


83 

84 


6 KOUNTap 

1 SLPMSatsNLNOi^-exLNQMJ 7(is»!LATM"EXUATM} 
Sl,P3NetSNLN(SN-BXl.NGN)/<EMUATN”6XL,AtNj 
C^3s(6NtNGM$SLPN3<)SNl,ATMI 
CSNs<eNLNGN?SLP5M»EJJLATN) 

XINTT<C3N-Cfj9}/{SLP,VSBSUPSN) 

yiNfi:{SURSN3,CMS»SLPNSt»CSN?/{S|,f?S[^,SLP,NS) 

NOUATsp 
K0UNT=K0UNT!1 
ir{KQUNT,rAT;iO> 30 TO 99 
DO 100 

100 I' f nNT,l,E,ARRAY(3i n ,AN0,YI^1T.GE,ARRAY(3, J*in .QO TO 101 
IPUOUNTfcGTV'lOl 50 TO 99 

NOLAT aj, 

GO TO 19? 

101 CAUU YriND(riMT,ARPAY(3, P , ARRAY t3, ! *fjl >» ARRAY ( g, U .ARRAY < 2 f I *1 > i 
, X ''I S } 

199 DO 200 I?1,]SJ~1 

200 IP(nNT,Ue,ARRY(3;U,AND»YlNT,GE,ARR5S?3, I + IM QQ TO 301 
ir{KOUNT,GTPlD) <30 TO 99 

N0t,ATs2 
GO TO g9? 

201 CaU YF4NIJ(?INT,ARRY<3. [J.ARRY<3«l>^l).ARRY<2ji),^^fjRY{2. K + l),XSN3 
299 IF(NOLAT,EQ:i} GO TO 305 

IPtNQLAT.eOJaJ 3C TO 310 
I^{ABS<XNS-^S^}^U, 0,0001) GO TO 4 OO 
IP<XNS.i.T,X8Ni 30 TO 300 

emuatn«xsm 

EMUNGNi?YlNT 

ENIAT.M»X(^S 

E(tUNGM=YINT 


QCl TO 1 

300 E>:UATN?X8N 
exUN6N=YlNT 
EXl-A,’TMfXNS 
£XI.NGM*YINT 
GO TQ 1 

305 IF{XShi.Gg.XTNT) Go TO 3Q7 
EMUATNcXSN 
ENL.NGN=Y1NT 
GO TO 1 

307 EXLATNtXSN 
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85 EyUN(SN»YlNT 

86 flO TO 1 

S? 310 ir<XNS,Q5,X!THT) GO TO- 311 ‘ 

88 eyiATMsXNS 

8? EX1-NGM^:YINT 

90 GO TO 1 

91 311 E^!LATMsyNS 

92 EMLNGM^YINT 

93 GO TO 1 

95 500 XCROSS*UNS?5?SN?/2;0 

95 YCROSSaYINT 

96 RETURN 

97 99 INOCpl 

9S return 

99 500 XcR0SS9RXI,AT(J 

?08 YCROSSxE;>{L,NGN 

3,01 return 

1Q2 501 XCROSSagXUA.TH 

?03 YcROSS*|=XUNQM 

g ^05 RETURN 

3,05 end 

1470 equality OR NON-EQU AL I TY COMPARISON MAY NqT BE MEANINGFUL IN LOGICAL IF EXPRESSIONS 



^ W ^ 
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SOtJRQyTlNH 1{f I NdUVTWOI, T^Oa, three UTHRIfEIi/THRES) 
DX»a-.TW(3tU{rMOU-TW3J) 

DY#TyREEtl»TflftSEl 
AOOTOYsDX^tD? 

THRESaTHftSgl^ADDTOY 
return 
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1 

2 

% 


to 


4 

y 

6 

7 

8 

? 

10 

11 

12 
1? 
14 
1 ® 
16 

17 

18 
1 ? 
20 
21 
22 

n 

2? 

2? 

28 

2? 

3? 

3) 

3? 

3| 


34 

37 




CGHTIM 

_SURRQUTINe SPinM(XcR03S,yoR0SSiGi?QJATUTj) 

C0MMQN/AI,T/St , array M f icon), NJ^Af?RYi^*1000 ) 

IF(ARRAY(2»8iU ;eo,XCRQSS^ Gq TO !?7 
IFURRAY<2',i) ,EQ,XCR0SS) 00 -TO 56 
NE p NNl 
DO 4$ n?l,f3E 

. lF''C(ARRAX<2;:iN)-XnR0SS)«{ARR'AVC^.frINtS}«XCR0SS} ,LT»0.) 00 To 55 
49 CONTINUE 
55 ApXCR0SS»»AR6aY{2i IN? 

BsARRAX(2#IS<’1)"ArrAY(2,IN? 

C*ARRAY(4» IS + 1)-ARRAY(4, IfO 
OpARRAni# IN + D-ARRAYU, IN) 

Gl»nAftC)/B}+ARRAY{4, IN) 

Tl«UA4DI/Bl + ARPAmi IN) 

NA * NJ”! 

DO &3 INfliBA 

IFC<ARRY<2,tN)-XCROSS)f»(ARHY(2iIN + l)i:XCRQSS);LT,0. ) 60 TO 42 
53 continue 
42 As5<CG05S«ARBY(2, IN) 

B = ARRY<2nNH>iARRY?2» IN? 

C = ARRY(4J:IN»1)1aRRYC4iIN) 

DrARRY(U1N?1)-ARRYC1»IN) 

GJ« (A<»C/B)?ARRY(4,‘|N) 

TU'*nA?D)/BI+ARRY(lin-) 

return 

57 GI«ARRAY(4,pn> 

TI“ARRAYU,R)I) 

QJSQI 

TJ«TI 

RETURN 

58 6I«AfilRAY(4,l) 

TPARRAYdilJ 

GJaGI 

TJ“TI 

RETURN 

END f 

1470 equality or NON-fiOUALlTY COMPARISON tIfiY NOT. BE MEANINGFUL IN LOGICAL IF EXPRESSIONS 





OODQOt){IODDSOOO(IOOGDO~DDOSDI10000000tltlBHOHII90ISOOO(IOIIOOIIOOaOOI)DOOaOOGOODODII9(IOOO 

I } ) 4 ] t II I II 111 ! II 14 IS II nil II II II ii n II » » II n II II II mi » IS » ii ii ii <i «i tt is h is it it ii ii si si s> si si ss si si ii » ti ii ii ii ii is it ii u is ii ii ii ii m is it ii ii n ii 

I I It 1 1 1 11 1 1 1 n 11 n 1 1 1 1 1 1 n 1 ti n 1 1 11 11 It 11 1 1 n ] n 1 11 1 M n 1 n 1 1 1 M II n II n 1 1 1 n 11 

2222222222~2222222""2222“222222222m22222t2222222222222222222B22222222222272222 
Ts 3 3 3 3 3 ”3 33 3 3 3 3 ~3 3 3 3 ~~ 3 3 3 3 3 3 3 3 3 3 3 3 3 3 n 3 3 3 3 3 3 3 a 3 2 3 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 3 3 3 3 3 3 3 3 a 3 2 3 

4 4444 4 4 4 4 4 4 4 4 M 4 4 4 U 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 M 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 14 444 4 4 4 4 U 4 4 4 4 4 4 4 4 4 4 4 4 4 4 44 

555SS$S55r5S5SS35~SS5S5~S55555555S55S5S55555SS5SSSSS5SS5S5S5555S55S55555S5S5S555 
E 6 6.6 6 S 6 G fi S G 3 { 6 8 S B G 8 6 G B 6 G 6 B G G S 6 6 S fi 6 6 S 3 G 6 6 B e G 6 S 6 G 6 S 6 S G G fi e G 6 E 6 G S S H G B S G 6 G B 6 G 6 e e G 6 e G 

m n m mn m 1 m m 1 m m m m n m j n n m 1 n m m m in m m m 1 m m 


99999399r999G9G9''99999~9993933993939399993999999G9399399399393S6G!iS9399399S99993 

I t ] i s t I I i lii II iMj u )) u ini a ;g it ii u n ?i n it a 3i i: ii it is it ii ii » tc 4i n a h ts «e v <i a sa si s? s si ss s( si it » ti i. a a ti (s ii ti u » ii n n ii u is k ii ii is to 

_____ MFcr-'Soai 



E-73^5 02 03-16-77 17.013 


1 

2 

3 

4 
9 
6 
7 

? 

i« 

11 

12 

13 

14 

15 

16 
17 

la 

1? 


22 

23 

24 
23 
26 
27 

29 
2? 

30 

31 

32 

33 

34 
33 

36 

37 
3B 
3? 

40 

41 

42 


CSOLVE 

. .subroutine BOLVen IADDR,§NREV,NSREV.,SUMi;sUM2,Deu1’A,B8US,KNSTR, 

, ARRAY) 

DIUBNSION A8RAYc4;U , ' 

DInENSlON S1C)MA(2) .fiBUS<UiDATA48)iSSHT<5) 

Q cohmqns south • north 

D0U9LE precision RUMl ( 1 ) i 8UH2 ( 1 ) »»DE;L^ A ( 1 } , ARCPAR C 6 ) » CQMPAK ( 3 ) 

inteuer SNTAPE 

integer ARcH 0 S( 6 ) ,GTAP,ScRB,SNREV(l)<,l!jSREV(lS , lADOR (1) »KNSTRU> 
INTEGER q,R 

integer REVtJd/DATAPE^DOUT 
integer 9CRC 
INTEGER DUHii.DUf^2;DUM3 
data DUHl/y/,DUf^2/0/# DUM3/0/ 

C0f1HQM/TAPEg/INTP;i0Ur,8TAPiNSTAPE^SNTAPE,SCRBiSCRC» DATAPSiDOUT 
COHHON/KQNTSi./NOPaRiNDNSiHOSNi IKNSTR,^gDlT,MAX 
EQUIVaUENCF: USN»DATA(1> ) , ( INSdUATa;?) ) 

DATA ARCN0S^1,2,3;3»0/ 

DA1*A SlQMA/2<>l.n/ 

data AVE/0.g/,SDEV/0.0/,KOyNT/0/ 

data intpos^o/ 
wRnEnogT>2i5oi) { iaddri pi 

2001 F0RMAT(////^^1XPSTARTJNG ADDR&SSeS 0P ARRAYS’i//'*3XPSNREVh5X, 
,0l2,/.3XPNBBEV'i'>X»O12,/,3Xl'SUMlM1X,012»/J’3XPSUM2P6X,012,/, 

,3X»«'OELTAt ,5X»012»/»3X| !8-3{AS’;3X>Q12,./,3 Xi ' KNSTR ' » 5X * 012 » 
,/,3xp aRRAY!;5X.012) 

RETURN 
entrx solve 

JTAPS19 

rehind sorb 

^J'<a^^^lKNSTB 

READ(SCRB) H, I2i (NSR = V( 1 U iniNONS ) * ( SNReV ( I ) i Pi , NOSN) i 
, (KNSTR( n , Pi, I jK) 

NDlM?NOPARtttNOSN + U 

CALL CUEAR^BOUS^NOIM) 

call FSTlHl{SUMliSUfl2,DELTA,BBlA§,ARRN0S,ARCPARfM3PAR.N0SN) 
AROPARU)=i:0DQ 
CO'HPAR<iJs-1.odo 
100 reao(scrb,end=2oO) data 
irnSN.LT.Ol GO To 150 
iarc=ins 

AR3PAR(2JsDATAr/> 
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43 


ARCpAR{3j8DftTA(7)»fr2 

44 


C0MPAR(2)a-DATA(8) 

49 


C0HPAR<3)s*DaTA(8)1h*2 

4^ 


RESI0sDATA<?)-^DaTa(6) 

4? 


ir<A8S(R8SlDJ .LE'iEDIT) 

4^ 


,CAUU eST^M{^j9PAR^^ISN,SOMPAR,S]:SMA^RE?ID) 

4V 


QO TO 100 

sa 

150 

ICONtO 

51 


IJKsNUM 9ER(SSREVUARC),KNSTR, iKNgTR) 

52 


IF( NK,QT.0J IcaN?KNSTR(): JK + IKNSTR) 

53 


CALL ESTfMAHCON) 

54 


CAUL ClEAR<B0IAS,MOPAR) 

55 


GO TO 100 

56 

200 

CAUU 6STIMC>SNREV;kNSTR, IKNSTR) 

57 


DO 500 K9N=1TN0SN 

68 


REWIND SORB 

5? 


READ<SdRB) Tl 

60 


IBCbNOBAR+1 

61 


IF{K8N,NBa) QO To 501 

62 


WRiTEMOUTiU) 

63 


WRIT6{6TAPiS5> 



WRlTEl^JTAPiJS) 

6i 


WRlTfi(GTAPil61 

66 


WR1T6(JTAP,16) 

67 


WRnS(GTAPil7> 

68 


WRiTBliJTAP.gO), 

6? 


WRlTe(QTAP'iJ9> 

78 


WR1TE(JTAP,19) 

71 


REWIND NSTAge 

78 

501 

CONTINye 

73 


DO 4Q0 IARCf!i,N0N3 

74 


CALL UPDATE 

75 


IF<KSN..NB,li QO TO 502 

76 


IREVtNSRBV(!IARC)/100000 

77 


ISEG?(NSREVllARC>-IREV« 100900 )<i:ig 

79 


WRITE(IOUT,1010) IREV,ISEG, (BBlASCPelPl'.NOPAR) 

7? 

205 

READ(NSTAPE} REVNOiNI, ( (ARRAY(C;R) ,Cfi$,4) ,r=i,n 

80 


lF(REVNO,NESNSReV(lAR:) ) GO TO 205 

81 


WRIT6(D0UT) DUMl'iDUR2*NIi lSEG,iRgV 

82 


DC ?07 KtIiNI 

83 


array (1,K)«?ARRAY(1,K)»60000, 

84 


ARRAY(2,K)pARRAY(2«K)»3600000, 
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85 


o 

c^ 


87 

80 

8 ? 

90 

91 

92 

93 

94 
93 

96 

97 
90 


9? 

100 

101 
102 
103 
1CI4 
103 
106 
107 
1U3 
10 ? 
110 


111 

112 

113 

114 
U3 
116 
117 


119 
11 ? 
'12 0 
121 


122 

12? 

121 

123 

126 


AFRAID (3, K)=l\8RAY{3it^) *3600000 , 

ITiMExARRAYh.K J 
ILATrARRAY(2«iK) 

IL0N=ARRAY(3;K) 

GO TO <2P8,S09i'2l0) iNDPAR 

208 ARRAYH,K>aiSRRAYM,K)=BBlAS(l> 

GO TO 206 

209 array (4,K) = liRRAY t4',(<)«{BBIAS(2)ftARRAlHtK) *BBIASU) ) 

GO TO 206 

210 XB8IASxR^IAg(3)^^ARRAYU,K )tttt2i+BBlASi2)oARRAYU,K )+BBIaSU) 
ARRAY(4,K)=ftRRAYU,iO«X8BlAS 

206 ARRAY (4,KJ=ARRAY (4, K) ttlOOOOOO , 

IFEAHTxaRRAXHiKI 

207 HplT£(DoUT> niME/IUATilUQNi lSE:AdT»DpM3 
ENDFILE DOUJ 

502 continue 

300 RFAD(SCRa,E5|Dx41C,) DATA 
ir( IS n.RQ,-!) go TO 400 
U (KSN.Nb,i;ANn. lSN,Nc,K8N> GO TO 3Qp 
I2=lBc+MqpA8«nsn-i) 

SSHTdlfiDATAO) 

SGHT(?)3pATA(6) 

DO 310 1*1,ROPaR 

SnHTCl)«SSHTa)-RBIAS( U<^DATA{7)«^^( Ij:l > 
SGHTt2)»SSHT(2)“FlBlASU2y«DATA<S)<»^^UTl) 

I?=I2+1 
310 CONTINUE 

SGHT(3)a<SSHT(l?'*'SSHT<2) 5 /2, 0 
SGHT{4>3SSHT(H-RSHT(3) 

SGHT(5>»sSHT{2)-‘oSHT(3) 

ir(A8S(0ATAt5)-DATA(6>) ;aT,EDIT) GO J,0 330 
A\>ExAVE + SSHT(4) 

SOEVtSDEV + SShT{ 4)fl^>2 
KOUNTxKOUNT*! 

330 irUSN.EU.lJ WR iTfrtGTAP, 1000 }NSREV(INS)/l000Q0,SfiREVUSN)7J 00000,.; 
,DATA(3),DATft(4),RBl ASa).BBlAS(l2^NOeAR)'i‘DATA(5) ,DATA(6). 

. (SSHT(K) ».K = i;'4) 

irnSN.FO,K3N)WRlTE< JTAP, 1000 ) SNREV(1,SN} 7100000, NSREV( INS) /lOOQOi), 
,DATA(3),DATA(4) ,BBI AS(l2>^N0pAR);BBlAgU)iDATA{6),DATA(5),SSHT<2> . 
,S^HT(1),SSHT(3) ,SSHT(5) 

QO TO 300 

400 continue 
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kO 

^-1 


^1 

og 



127 

129 
12 ? 

130 

131 

132 

133 

134 

135 

136 

137 

139 
13? 

140 

141 

142 

143 

144 

145 

146 

147 

149 
14? 

150 

151 

152 

153 

154 

155 

156 

157 
153 
15? 
160 
161 
162 

163 

164 

165 

166 

167 

168 


410 coNTiNye 
500 RcWlND SPRa 

WRlTEaoUT.lO) 

RSWIND SNTAf!5 

DO 420 iniBOSN 

IEsIBC*fviQpA54l 

nEV¥SNRRVa)/lO00OQ 

I3EGt(S.MREVM>-IREV<H0OOOO>/1O . „ 

WRITBUOgTiiOlO) IREV, IStGi (BB IE) 

405 RrAD<SNTAPE5 RE VNO i NJ , ( ( ARRAY ( Ci;R ) , Cfil f 4 ) , Rsi, NJ ) 
ir<REVNOfMe;s,MREV(l) ) GO TO 40? 

WRITE(DOUT) DUMl.DUM2»NJnSEG, IREV 
DO 407 K?1 iSJ 

ARRAYUiK)«ABRAYa,K><^60000, 

ARRAy(2.K>sARRAY(2iK)<*3600 0 00, 

ARRAY(3iK)=ABRAY(3iK)»360OQO0. 

ITlME = ARRAYh»KJ 
ILAT?ARRAY(2;K) 

IL0N?ARRAY{3;K) 

GO TO C408,509,4U)iNOPAR 
4 08 ARRAY<4,K)=ARRAY(4iK>e08IASa0C) 

GO TO 406 

40 9 ARRAy<4,K)sAgRAY(4,K)«<03IASnBC^!l)»ARRAyU>K)+'30lAS< IBC)? 

GO TO 406 

411 YBlASaQ9lASUaC + 2)ttARRAY(UK7^M>2^ + BBHASnBQ*l><»ARRAYU»K) 
Y8 IASsXBIAS»B 81 AS( tBC ) 

ARRAY(4,)^)=fAgRAY{4»iO «YBIAS 
4 06 ARRAY<4,K) = ARRAY(4»><)»1000000, 

ISEAUTf ARRAn4,K) 

407 WRITe(DQUT) I Tj ^E I L AT » lUON , J SB AyT ^DOMS 
ENDFIUE 0QUT 
IBCsIBD*NOPAR 
420 CONTINUE 

ave=avb/koust . 

SDEVirSDRV-AVB^vtt? 

SDEVtS9RT{SD6V/K0UMT) 

WRITE(GTAP,«8) AVE 

28 FORMAT(«0»,iyEAN OP RES I DUALS U5X'; F1905 ) 

KRITE(GTAP,29) SOEV 

29 FORMAT(«0« iISTArDARD DEV I AT ION ' » , FlO, 5 ) 

11 FQRMAT< flJ,l7X, 16X, 1N0RTH"TQ*S0UTH REVS*,/ 
,»REV NOi,5X7,iSEG NO i * IIX , » COEFF 1 0 1 EN J.S i ) 
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16? 

170 

171 


172 

173 

174 


17? 

176 

177 


176 

17? 


18Q 


181 

182 


VO 

oo 


183 

184 
18? 
186 
187 

tt ^ tt W 


15 rORMATdl* ,63Xi » A PR I DR I M3X * ' AO JUSIBD < ) 

16 F0RMAT<lX, iSEVOLUTIONSSdX, 'CROSSING RTS'. > ,12X, •BtASES* ,9X, 

,»SEA SUrKACB HE1QHTSS3XPSEA SURFACB HE | GHTS ' 1 3X AVERAGE », 3X » 

, 'RESIDUAUS'J 

17 r0RMAT(2X, »S/Sf,3XV'S/Nt ^4X, >LATITUDB LONG I TUOE » ,'6X , »N7S»,7X, 

, 'N/S ',7Xi I S^N » ;?X ; 1 N/3 8Xi; t S/N ' k 7X, *M/S1,5X;< (ADJUSTED) h4X» 
,'N/S') 

13 FORMATC'O' ,19X, ’SQUTHaTCsNORTH ReVSifc/,'REV N0»|5X,'SEG NO»,UX, 
, 'COEEFICl'ENTS' } 

19 f QRMAT(l4Xi!tDHGRE6S) t ,2X, ’ (QE6REES)**AXi 1 (METEPS ) ♦ , 2X, 

, ’ (METERS) »,3X, » (METERS) 'ISX, ' (MET.EKS ) i , 3X, ' (METERS )'# 2X, 

. ' (METERS)' »3X# ♦( METERS) * >3X, ((METERS n) 

20 F0RMAT<2Xi 'B/N' ,3X, 'N/S' IaX, 'LAT iTUDB LONG I TUDL ' , 6X , f S/N ' i 7X , 
PS7N'*7 x;'N^S'iVX;(S/v|*»8Xj » N/S ' »<7X r iS^N * » 5X>‘ ' (ADJUSTED )♦ i 4X » 

, »S7N') 

1000 F0RMAT1/12(1X,I5'> ,lXi2(lX,FlO,6)«2X,?nX»'F9,4) ,3UX/2(1 XiF9,4)>> 
1010 F0RMAT(( 'iH»7X, i4,1DX,E17,7»?X7e17^7,5Xi617»7) 

STOP 

END 

MEMORY' expanded, USE SLIMITS OR COREa OPTION FOR NEXT RUN 
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L8L CHARLYi5f5AHT SUBROUTINE CHARLY 

TTL SUBROUTINE QHARIY - DYNAMI5 MEMORY ALLOCATOR 

start OE documentation comment cards 

PROGRAM title - SUBROUTINE CHARLY 
PROf^RAM NOi 1X2293 

written by Dan chin (Wolf research and development) 

REWRITTEN BY TOM HARMON <NASA WALLOPS-AMS) 1/74 
MOOT, PIED BY TOM NORTHAM JWRPC) 5/74, 6/74 

HAROHARfel/SOFTWARE SUMMARY 

COMKOTER REQUIRED ** HONEYWELL 600/6000 SERIES 
SYS,TeM EXECUTIVE - GECOS m 
PROGRAM language - GMAP 
PERI.RHERALS - NONE 

PURROSE 

THIB SUBROUTINE DYNAMICALLY ALLOCATES ARRAY STORAGE FOR USE 
WITI|1 THE SEAHT program, 

HEiyOD 


WORD 31 (DECIMAL) IN THE SLAVE PROGRAM PREFIX IS SCANNED 
TO HND the location op available memory FOR ARRAY AULOCATI 
IP sufficient memory 1SN*T AVAILABLE! THEN AN ATTEMPT TO 
GEMpRE MORE MEMORY IS MADE, IF GCOS DENIES THE REQUEST, 
THEN CHARlY returns TO THE NON STANdARQ ERROR EXIT, 
otherwise, the array locations ARE STORED INTO THE SOLVEI 

calling sequence and a caul to solve! is executed, 
calling sequence 

call CHARLY, {KORBBi IMOREfNcOREj IC0RE,K0REG,'*NNN) 

KOREB » the amount OP CORE THAT MUST BE AVAILABLE FOR 
BUFFERS AFTER' CHARLY HAS BEEN CALLED, 

IMORg s THE SUM OF ARRAY SIZES, 

NC09E 8 NUMBER OF ARRAYS FOR WHICH CORE MUST BE 
ALLOCATED. 

■IPRE « number of words QF CORE REQUIRED BY .EACH ARRAY, 


charlyoi 

CHARLY02 

tt#CHARLY03 

«CHARLYo4 

ttCHARLYoB 

OCHARLY06 

^»CHaRLYo7 

ttCHARuYoa 

ttCHARLY09 

ttCHARLYlO 

ttCHARLYll 

ttCHARLYl2 

!>0HARLY13 

„CHARLYi4 

ttCHARLYlS 

«CHARLY16 

ttCHARLYl? 

ttCHARLYlB 

^^CHARLYi9 

<tCHARLY20 

ttCHARl,Y21 

<hCHARLY22 

«CHARLY23 

ttCHARLY24 

«CHARLY25 

0NCHARLY26 

ttCHARLY27 

ttCHARLY28 

#CHARLY29 

«CHARLY30 

«CHARLY31 

^^CHARLY32 

«CHARLY33 

tt0HARLY34 

«CHARLY35 

ttCHARLY36 

ttCHARLY37 

ttCHARLY39 

h»CHARLY39 

■ttCHARLY4Q 

<»CHARLY41 

OCHARLY42 



101 


17 . 01 ’; 


susROuriNe 

charuy n d5(hamic memory allocator Page 

3 

43 

0 

KOREG S THE AMOUNT OK CORE LEFT IN USER JOB SEGMENT 

*CHARLY43 

44 


IF POSITIVE, "THE AMOUNT OF CORE TO BE QEMORE-D 

■«CHARLY44 

45 


IF NEGATIVE, ' ' ' 

»CHARLY45 

46 


$NNN s RETURN WHEN QQR6 ALLOCATION NOT SUCESSFUL 

ftCHARLY46 

47 


■ 

ftCHARLY47 

43 

# 


*CHARLY48 

49 

ft' 

RETURNS . 

ftCHARLY49 

50 

IV. 

1 

ftCHARLYSO _ 

51 

ft 

normal Only '(TERMINATION OCCURS IN cONSOL 

ftCHARLY51 

52 

ft 

. ON 'INSUFFICIENT MEMORY AVAILABLE) 

#CHARLY52 

53 

ft 


«CHARLY53 

34 

ft 

INPUT 

«CHARLY54 

55 

ft 


ttCHARLYSS 

56 

ft 

VIA CALLING SEQUENCE 

*CHARLY56 

57 

ft 

• 

ftCHARLY57 

58 

ft 

OUTPUT 

ttCHARLY58 

59 

ft 


ftCHARLY59 

60 


chably Generates a simulated call sclvemargi,arg2,...) 

ftCHARLY60 

61 

ft 

WHEBE ARGl POINTS TO THE ADDRESS OF ARG2, AND ARG2 THRU ARQNftCHARLY61 

62 

ft 

CONTAIN the addresses OF THE DYNAMICALLY ALLOCATED ARRAYS, 

ftOHARLY62 . 

63 

ft 


ftCHARLY63 

64 

ft 

SUBRR0G8AMS REQUIRED 

«CHARLY64 

65 

ft 


#CHARLY65 

66 

ft 

SOLVE I 

ftCHARuY66 

67 

ft 


*CHARLY67 

68 

ft 

called by 

ftCHARLYsa. 

69 

ft 


«CHARLY69 

70 

ft 

seamt main program 

*CHARLY70 

71 



ftCHARLY71 

72 

ft 

external references 

ft(;HARLY72 

73 

ft 


ftCHARLY73 

74 

ft 

NONB 

. «CHARLY74 

75 

ft 

BIBLIOGSAPHY 

ftOHARLY75 

76 

ft 


*CHARLY76 

77 

ft 


#CHARLY77 

78 

ft 


#CHARLY78 

79 

ft 


«CHARLY79 

00 

<i,<>ftftft* end of documentation comment CAWdJ* «*ftftftftftft*ftft»ftftftftftft*ftftftft«ftft»ft*ftCHARLT8Q 



6 - 7 V 17, 017 subroutine CHARLY - DJ(NAMIC MEMORY 



QOOQQ^ 

82 

CHARLY 

SAVE 


000002710000 
00D?A76i0000 
000217754000 
0O0P17741OO0 
000002 7200 31 

010 

01(1 

010 

dm 

don 

63 

LXLu 

St 

000044 

7400 

00 

Gif) 

84 


STXQ 

X 

ooono? 

5500 

00 

DID 

85 


SBAR 

9*1 

000000 

2350 

07 

oon 

86 


LDA 

tttt,DL 

000777 

3750 

07 

oon 

87 


ANA 

A0777,DL 

ooooli 

7358 

00 

Oqii 

88 


ALS 

9 

0 1 0 0 0 1 

1750 

07 

oon 

89 


SBA 

liDL 

000037 

7510 

07 

oon 

90 


stca 

51*07 

4il» 

OOOQQ4 

2350 

33 

Ooo 

91 


LDA 

000006 

7350 

00 

000 

92 


ALS 

6 

000153 

7510 

06 

01') 

93 


STCA 

0TTAlY+1*D6 

000153 

2350 

00 

oin 

94 


LDA 

OTTALY*! 

000152 

7550 

00 

010 

95 


STA 

OTTAlY 

000005 

6200 

3l 

oon 

96 


EAXa 


OOh 51 

7400 

00 

Oin 

97 


STXO 

INTAlY 

Y o^o 
00 37 

CHaRLX02 
2350 00 

doij 

98 


LDA 

31 

777777 

3/50 

03 

oon 

99 


ANA 

-^1*DU 

OOOOOl 

0350 

03 

Don 

100 

AGAIN 

ADLA 

1*DU 

77777A 

3750 

03 

Don 

101 


ANA 

40777776*nU 

00(1152 

7550 

56 

oin 

102 


STA 

OTTALY, ID 

d00002 

b070 

04 

oon 

103 


TTF 

2, IC 

000036 

7100 

00 

Din 

104 


TRA 

CKMEM 

Oooisi 

7200 

96 

Oin 

105 


LXUU 

INTALY* ID 

000334 

74(jo 

00 

Din 

106 


STXO 

tt + 1 

000300 

0350 

03 

(ioo 

107 


ADLA 

*h*,DU 

000 J25 

7100 

00 

OlO 

lofl 


TRA 

again 

000151 

7200 

56 

010 

109 

CKMEM 

LXLO 

INTAlY, id 

000340 

7400 

00 

010 

110 


STXO 

fl + 1 

000300 

0350 

03 

doj 

111 


APLAi 

««,DU 

000337 

7510 

70 

000 

112 


STCA 

3li7o 

000337 

7200 

00 

Ooo 

113 


LXLO 

31 

030337 

1600 

0 0 

Don 

114 


SBXU 

31 

000300 

1600 

03 

001 

115 

•< 

SBXQ 


0003QO 

6350 

10 

ooo 

116 


EAA 

0,0 

Y CARD 
000022 

CHARLY03 
7310 00 

003 

117 


ARS 

18 


ALLOCATOR 


PAGE 4 

CHARLY82 


ACCESS KOREP 

STORE IN UPPER half OF INSTR AT LOG X 

' GET BAR IM AL 
NO, OF WORD BLOCKS 
MULTIPLY BY 512 
UPPER CORE ADDRESS 

NO. OF elements 


initialize working tally word 

ARRAY ADDRESS 


CHARLya3 
CHARLY84 
_ CHARLY85 
CHARLY86 
CHARLY87 
CHARLY88 
CHARLY89 
CHARLY90 
CHARLY91 
CHARLY92 
CHARLY93 
CHARLY94 
CHARLY95 
CHARLY96 
CHARLY97 


GET LOWER MEMORY ADDRESS IN AU 
MAKE ADDRESS EVEN 

and store Into callin'? seovence 

CHECK FOR TALLY RUNOUT 
YES 

SIZE OF ARRAY 


SIZE OF LAST array 



UPDATE LOWER LlMjT 
UPPER MEMORY LIMIT 
LOWER MEMORY LIMIT 

GET XO IN A'J 

SHIFT AU To AU W SIGN SET 








CHARLY98 

CHARLY99 

charlyoo 

CHARLYOI 

CHARLY02 

CHARLY03 

CHARLY04 

CHARLY05 

CHARLY 06 

CHARLY07 

CHARLY08 

CHARLY09 

CHARLYlO 

charlyu 

CHARLY12 

CHARLY13 

CHARLY14 

CHARLY19 

CHARLY16 


CHARLY17 
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00006 7550 31 

000 

iia 

STA 

6»1<5 

CHARLY18 

50LS4 6050 OQ 

01) 

ii« 

TPL 

PCALL WE HAVE ENOUGH MEMORY 

CHARLY19 

OOOOO 5310 00 

009 

120 3M0RE 

NEQ” ' 

MAKE DIFFERENCE POSITIVE 

CHARLY20 

U777 0350 07 

dOQ 

121 

ADLA 

1023, DL ROUND TO NEAREST iK 

CHARLY21 

00n2 7710 00 

003 

122 

ARL 

fo ■ ' . 

CHARLY22 

D0150 7550 00 

019 

123 

STA 

HEMSIZ 

CHARLY23 

D0061 7510 07 

019 

124 

STC« 

MME*1i 07 .STORE K REQi) IN MME ZERO WORD 

CHARLY24 

50012 7350 00 

doo 

125 

AL3 ... 

10 . 

CHARLY25_ 

50037 0550 00 

d09 

126 

ASA 

31 update upper memory limit 

CHARLY26 . 

50011 0010 00 

009 

127 MME 

MME: 

GEMORE ASK FOR MORS MEMORY 

0HARLY27 

50000 000000 

doo 

128 

ZERO 

d, tttt 

CHARLY28 

50'502 7100 04 

doo 

129 

TRA 

2.IC DENIAL 

CHARLY29 

50135 7100 00 

Old 

130 

TRA 

PCALLl 

CHARLY30 

50065 7l00 00 

dio 

131 denial 

TRA 

SMESqi 

CHARLY31 

50900701000 

030 

132 3MESG1 

CALL 

eONSOLlMEMSiZi IREPLY)DENRET 

CHARLY32 

50973710000 

013 





50217000204 

on 





50150000000 

Old 





CARD CHARLX04 






50'!-* 000000 

on 





)0JS 710000 

Oio 





)01 2340 00 

on 

133 

SZN 

ireply 

CHARLY33 

5(5131 6000 00 

010 

134 

TZE 

DENRET 

CHARLY34 

50143 0540 00 

on 

135 

APS 

GMRINO 

CHARUY35 

L0000701000 

030 

136 

CALL 

CLKTiMITIMZRO) 

0HARLY36 

50102710000 

010 





50217000210 

010 





50145000000 

dio 





30130 2350 00 

Old 

137 

LDA 

DENALI 

CHARLY37 

50064 7550 00 

Old 

138 

STA 

denial 

CHARLY38 

50060 7100 00 

Old 

189 

TRA 

MME 

CHARLY39 

10000701000 

03) 

140 3HeSG2 

CAUL 

CLKTlM(TIME) 

CHARLY40 

50111710000 

Old 





50217000214 

Old 





50146000000 

Old 





50146 2350 00 

019 

141 

LDA 

TIME 

CHARLY41 

5Q145 1750 00 

Old 

142 

SBA 

timzro 

GHARLY42 

5GD36 1150 07 

000 

143 

CMPA 

30, DL 

CHARLY43 

CARD CHARL>:05 






50560 7100 00 

010 

144 

TRA 

MME 

CHARLY44 

30000701000 

030 

145 

call 

EONMeS 

CHARLY45 

50120710000 

Old 

■* 
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. SUbROUTINS CHARLY , D^(MAMIC MfiMORY 


00?17Q0Q221 OIQ 

10000701000 033 

00124710000 010 

00217000222 010 

00143000000 010 

00144 0540 00 OlO 
00144 2340 00 010 

00131 6000 00 dio 
d0060 7100 00 010 

00105 7100 UD QID 
00013i 

00717221051 on 

00001161003 000 

00217741000 QIC 

00001710000 QIO 


.1^6 -CAUL. . eiKTiMnriZRg 


147 APS 

148 SZN 

149. . rZfc . 

150 TRA 

151 PENAl,! TRA 

152 PENRei RETURN 


C;,KCNT 

CLKCNT 

DENRET_ 

MME 

SMESG2 

PHARLY,! 


03143 2340 OO 

dip 

153 

fCALLl 

SZN 

GMRIND 

OOL54 6000 00 

013 

154 


TZE 

FCALL 

C<^RD CHARt.X06 






20 ‘" 0701003 

030 

155 


CALL 

CONGRT 

00 S 2710000 

OIQ 





00 '000231 

013 


. 



00154 7100 00 

dip 

156 


TRA 

FCALL 

00300000000 

QOO 

157 

GMRIND 

Df;C 

d 

77777777766 

000 

153 

clkcnt 

DEC 

-10 

00014? 


159 

T1M2R0 

bss 

i 

000146 


160 

TIME 

RSS 

1 

030141 


161 

IREPLY 

3SS 

1 

000150 


162 

'•^EMSIZ 

BSS 

f 

00000 OOOO 00 

doo'^ 

163 

INTALY 

tally 


00000 0000 00 

OOP 

164 

OTTALY 

TALLY 

(M> , 

00160 OOOg 00 

Old 

165 

} 

tally 

FCALL* 



166 


symref 

SOLVEl 

50000 7010 00 

030 

167 

FCALL 

TSXl 

SOLVE! 

00216 7100 UO 

dio 

168 


TRA 

mcall 

00217 000000 

010 

169 


ZERO 

?EiL. . 

00160 0000 00 

013 

170 


ARQ 

»*1 

OOQ160 


171 


BSS 

30 

00001710000 

010 

172 

VCALL 

RETURN 

CHARLY 


AG'^ 


00000000000 


000 


ALLOCATOR 

t 


PaQE 6 


_.CHARLY46_. 


CHARLY47 

CHARLY48 

CHARLY49. 

CHARLY50 

CHARLY51 

CHARLY52 




INPUT tally word 
WORKING TALLY WORD 


MAXIMUM 0” 30 ARRAYS 


CHARLY53 

CHARLY54 

CHARLY55 


CHARLY56 

CHARLY57 

CHARLY58 

CHARLY59 

CHARLY60 

CHARLY61 

CHARLY62 

CHARLY63 

CHARLY64 

0HARLY65 

CHARLY66 . 

CHARLY67 

CHARLY68 

CHARLY69 

CHARLY70 

CHARLY71 

CHARLY72 . 



105 


SUBROUTING PHARUY ^ PYNAMIC MEMORY ALLOCATOR 
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000220 233Q2l5lfl370 000 

GND of BlfJARY CARD CHARL)^07 

173 gNp 

22? IS TfHE next available LdCATlON, 

JMPA 730601/052373 UMPR 730601/092373 

THeRI WERE MO WARSING FLAGS IN THE ABOVE ASSEMBLY 


JMPC 730601/052373 
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1 


4 

J 

6 

7 


13 

« « « V W 
11 
12 
ft fr jf i«l 

1<* 

ft ft ft ^ ft w 

13 

16 

M ft^ftW 

ON 1 7 

ft ft ft ^ ft W 

le 

1? 

ftftftf^ft w 
28 

21 

ft'^ft^ftN 


22 

23 

24 

23 

26 

27 

2§ 

2 ? 

30 

31 

32 

33 

34 

35 


CCOMSOL 

BONSOU.^MEMSZEnREPLyi^*) 

dimension MES(lO),M£S0ai)iMESlfi;iyit1fiS2<8)'iMES2l<9)iMGS3(9h' 

, M5S3(f) 

CHARACTBRtt6 «ES»MESC,''^GS1,MES2*M&S21«HES3 ;mES5' 

CHARACTEH«6 ftONSLW;NA'JSlVNANS2#NANS3#NANS4 
CHARACTEB*6 NANS5,‘NANSt>, NANS7,NANS8 
, character *6 IZERO 
CHARACT£Rft6 NRcPLY 

data MES/60HPLEASe Set other ^OBS UR®c to oo so they will swap out 
4U2 data STATEMEMTi VARIABLE LIST EXCEEDS HTgRAL LIST 

I • • 1 f 

DATA MES0/6SHTVPE Y6S ft* A'ND SET ALU OTHER ftft URQC ft* TO 00 OR 
402 data STATEMbNTi triable LIST E^^CEEDS LITERAL LIST 
,#« NO ftft f 

DATa MES1/66HGEODYN NEEDS K ADDITIONAL CORE TOR A TOTaU OF 
4U2 data statement, VARIABLE LIST EXCEEDS LITERAL LIST 


• K » • I / 

data MES2/40HTYPE YES AND TEMPORARILY SWAP OTHER JOBS OUT,,, / 

402 data statement, VaRJABle LIST eXcEEDS LlTpRAL LIST 

data MBS21/54H BY 3ETTINQ ALL OTHER URGENCIES TO ZERO,,, OR,, 
402 data statement, VARIABLE LIST EXCEEDS LlTgRAL LIST 

*DATa MES3/54HTVPE NO AND GEOOYN WILL STOP EXECUTION aND ENdJOb, 

402 DATA STATEMENT, VARIABLE LIST EXCEEDS LITERAL LIST 

data MES§/5S?(ftftGeODYNfttt GOT THE GORE IT NEEDED,., THANK YOU, 

402 data statement, VARIABLE LIST EXcEEDS LITERAL LIST 
» / 

data CQNSUW^6H0O00T//,NANSl/6HyE8 Y.NANS2/6HY 
,NANS3/6HN0 /»NANS4/6HN / 

data NaNS5/6hYFSO 00/» ^ ANS6/6HYOOqOO/*NANS7/6HNOOOOCI/, NANSS/6HNOOOO 

, 0 / 


data ize«o/6Hoooooo/ 

NREPLY*I2ERQ 
CALL MEHSIZtMTSlZE) 

MTSIZErNTSIie+MEMSZE 
PRINT 1OO,MB0S7.E,MTSIZE 
100 FORMATaX,0l2,lX,0l2) 

IF(MTSIZE,GT,80.0R,,NE'1SZE,GT,25? GO ,T0 15 
ENCODE {MESJt3),10) MzHSZE 
ENCODE (MESJUn, ,1D) MTSIZE 


CONSOLOl 

C0NS0LQ2 

C0NS0L03 

C0NSOL04 

CONSOL05 

COnSoL06 

C0NS0LQ7 

C0N50LQ8 

CONSOL09 

CONSOLIO 

CONSOLll 

C0NS0L12 

C0NS0L13 

C0NS0L14 

C0NS0L15 

C0NS0L16 

CONSOLl? 

CONSOL18 

C0NS0L19 

C0WS0L2D 

C0N50L21 

CONSOL 22 
C0NS0L23 
C0NS0L24 
C0NSOU25 
C0NSQL26 
C0NSOL27 
C0NS0L28 
C0NS0L29 
CONSOL30 
C0NS0L31 
C0NS0L32 
C0NS0L33 
C0NS0L34 
C0NSOL35 
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O 

--4 


36 

10 

37 


39 

a 

3? 


40 


41 


42 


43 


44 


45 


46 


4? 


49 

12 

4? 


50 

13 

51 

14 

52 

15 

53 


5 ^ 


55 


56 


57 


59 


5? 



F0RMATU6) 

CALL CNSLlO!tJOMSLWi‘HE31|U> 

CALL’ CNSLI0»C0MSLW(ME32,8) 

CALL CNSLI01C0MSLW,'ME321.,9) 

CALL CNSLIOICONSLW,MES3,9) 

CALL CNSLIOtCONSLW»ME30,U^NREl?LY) 

IF{NRERL55.E9iMANSl,OR,NREPLYiEQ4NANS3J GO TO 12 

IF(NREBi.<.Ea'^NANS5,0R,NREPLY-,EQ.f NANSfi} GO TO 12 

LF(NRERLY,ESiNAN'S3,OR,.'jREPLY.EQ*l!JANS8) GQ TO i? 

IF(NREBL)5.E®:NANS7,OR,NREPLY,EQ.NANSg} GO TO 13 

NREPLYriZERQ 

GO TO 11 

IREPLY=1 

GO TO 14 

IREPLY=0 

return 

return 1 

entry gonnes 

CALL CN3LIOtCONSLW|ME3,10) 

RETURN 

entry conqrt 

CALL CNSL10tCOMSLW,ME35,9) 

return 

END 


C0NSOL36 

C0N3OL37 

CONSOL38 

C0NS0L39 

C0NSOL40 

C0NSOL41 

C0NS0L42 

C0NS0L43 

CQNS0L44 

C0NSOU45 

C0NSOL46 

C0NS0L47 

CONSOL43 

CONSOL49 

CQNSOL50 

CONSOL51 

CONSOL52 

C0NS0L53 

C0NS0L54 

CONbOLSS 

C0NS0U56 

C0NSOL57 

CONbOUSa 

C0NS0L39 
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CCLKTIM 

SUBROUTiNli RLKTIM({TI.'1B) 
■ CAUL PTIHE(ftTIM6) 
ITIHE=ATIME«3600, 

return 

END 



QOOOOtlOOOOOOOQCOC090''OOtlOOOOSOOOQOBOODDDtlQinillOBDOI100llOOQSOODDDOtlOOOOQ9l!Ot)liDO(tin)0 

I 1 I 4 i ( II I II nil n 14 » II till nil ii iiiiii »»ii iiii uii »u» is min Xihkiuh isiiii ii <i si si sisisi ssusisi siiiii uiiii isisiiiiii laii ii n n isii ii ii iiii 

1 1 1 1 1 1 1 1 1 f I n 11 1 1 11 1 1 1 1 n 1 11 1 11 11 It II n n 1 1 1 1 1 1 1 n 11 M n n 1 1 M 1 1 n 1 11 M It n 1 1 1) 1 
2222222222“2222222“22m~22222222222222222222222222n22Z2222222222222222222n212 

'■333333"3333333"3333’"'333"3333333333S33333333333333333333333J3333333333J33333333 
4444444444444444444~444444444444<444444U144M4U444U44444444444444444444444441 
SSSS55S5S~5555355~5S555''5S55S55SS5S5S55SS5S535SSS55S555S5S555555555S5535S5555SS3 
B66S6e66&SES6B6GB668Gfi6CEBGSEBS6D66BC6E66G68E6SS8ES6B8686S66(HGE366B636ES6GGS66 
77771777717777777777777711177117777777117177177777771777771177777777777777777777 
8 8 S 8 t 8 8 8 8 E 8 8 B 8 8 B S 8 8 ~ 8 8 8 E 8 8 E 8 8 B a G B t B 8 8 B 8 8 8 B 8 8 8 S E 8 3 n 8 8 8 e 8 3 B G t 8 I ( 8 B 8 B E 8 8 a G 8 8 8 8 B a 8 
9993339E''9999939~99999~93999399999939999S939999999999993999999999li939999E3999399 

■ I I I s I ) I 1 II n 1)11 II IS II It II IS >1 II 11 )i )i IS II II II ;i :i <i » ii i« is is ii it is u <i » ti ii is ii ii ii is si si s) >i <i ss si si si si n i. it ii ti is ii ii si ei ii ii tt ii it is n ti it is ii 
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i 

4 

4 

4 

3 

6 

7 

? 

u 

12 

13 

11 

13 

16 

X7 

18 


21 

22 
23 
26 

25 

26 

27 

28 
2 ? 
38 

31 

32 

33 

34 
3^ 

36 

37 

38 
3V 
48 

41 

42 


CESTUI ; 

SUaKCUTINE 8STIMl{SUMl,SUM2,Dei,T/i,PBUSiPARN0liPARTLnNDXl,MDX2) 
poUaUE PRECISION SUMi(t>iSUM2(lUDEL,TSa),PARTL.I (1 ) , PARTLC a ? , 
Dimension bbia'^cd ^sigmauj ■' . 
integer RARROI a);PARN0C'i'SCRC,3NREVUJ ,KNSTR(-1) 

LOGICAU NEWA8C 

COM(4GN/TAPEg/IHTP;iOUr,GTAP,NSTAPEiS^TAPE,SCRBiSCRC 

data NRWARC^ 'TRUE./ 

data BIaS$G/'0,D01;.2i10,^' 

lNDXNOU>’NDlM*KNl)-CH»n"-l> >/2 

NPARl=NnXl 

MPC0MBNPAR14NDX2 

NDlMsNPARl+SPCOM 

nstart?nraRi+i 

IEND = NPARI«^2 
I1»INDXNQ(N0IM>+NdIM 
CALL CLEAR(guMl»2ftIl) 

CALL CLEAR(gUM?,2ttNDM> 

return 

entry ESTlMfcPAPNOciPARTLC|SlGMA»RESJD) 
lF(SIGMA(l)-?LE.q.0,OR,SIGMA(2).U&,0,P J RETURN 
1F( fNOT.NEWARCI GO TO 25 
NSTORE?INDX80(NPARn + ^DIM 
CALL CL5AR(8UMl,2ttNSTDRE) 

CALL CLEAR < SUMP. 2^^NPAR I > 

NEHARCf, false, 

25 DO 50 lal.NSARI 
IIsHNPARI 

PARTLI { U)=aARTLC(I ) 

50 PaRNOI ( ll)apARNOC+l 

WT=1,0/(S1gHa(J )<m>2 + SICMA<2)Im>2) 

DO 75 lal, IEND 
Il = PARNnIU 1 
T = PARTUHn«WT 
SUM2ni)fSUM2( ID+ToRcSID 
DO 65 Jaf , ISND 
ISTsINDXNO(U)+PARMOI (vJ) 

65 SUMlUST)sSUMinST)*T*PARTLlt J) 

75 CONTINUE 
RETURN 

entry ESTIMAdCON) 

IFdCON.EO.^) GO TO 3 OQ 
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43 

44 
4^ 

46 

47 
4$ 

4? 

5(J 

51 

52 

53 

54 

55 


5? 

5? 

66 

61 

6 ? 

63 

64 
6? 
66 
67 
6 8 
6? 
76 

71 

72 

73 

74 

75 

76 

77 
7§ 
7? 
60 
61 
02 

83 

84 


DO 85 1BIC0S,MPARI 
ji-iWDiJNa<n + i 

35 SUMUIi)¥5UMl(iij+l,0Dp78IASSG<]}»42 
_100 CALL Sy'1lNV{SUMl,ND}M^NPARI,DEUTA) 

u=ino^nq<nstart) ■■ * ' 

DO 200 ItNST.ART,NDIM 

u*0 

DO 160 UtIiSPARI 
... PELTA{L,)?0*^DO 
IST#0 
DO 140 
41»1ST*U 
vt2=lST^l 

DELTA(U>?DEliTAa>+SUMif,Jl)4SyMM02? 
140 ISTsIST4.ND1H«J 

IF(L,‘,Ea,NPARn GO TO 150 
' LPlsU+l 

00 145 4?LPS;NPAR1 
jZ^ISTtl 

0EUTAa)*DEli’TA(L)*SgHl(Jt}«SUMl{k^2) 
145 |STs1ST*NDIM4j 

ISO SUH2U>«S0M2{n-SUMa(Lr}ttDHUTAU,l 
DO 155 KfliSQIn 
KlsJl+K 
.K2sLl*K 

155 SUMl(Kl)=iiSUflt(KlJ-DEl,TAa>ttSyMl(K2) 
160 LlBU+NOiM-i: 

Jl=l 

DO 175 J?1i!3RAPI 
SOHlUDPDEClTAtJ) 

175 glsgi+NOlM-’iJ 
200 I1«U+ND1M-1 
iSTaO 

DO 250 I*1,SPARI 

DEtTA(I)?0,5D0 

I1»I 

DO 240 Jfl.l 

OELTA<l>*DEa'TA( J )+SUMi< 11?4SUM2{U3 
240 Jl = U + NDlM-8 

lp{UEQ,NPABn DO TO 248 
1P1=I+1 



Zl I 
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86 

87 

8? 

8? 

90 

91 

92 

93 

94 
9» 

96 

97 

98 
9? 

100 

101 

102 

10^ 

104 

105 

106 

107 

108 
10 ? 
110 
111 
na 

U3 

114 

lis 

116 

117 

110 

11 ? 

420 

121 

122 

123 

124 

125 

126 


DO 245 J¥IP1;nPARI 
li^ISTikJ 

245 DELTA( n?DEe;TA< 1 )+SUMl( Il)ttSUM3<vJ} 

246 B8lASn))?e8lAStI)+D6UrA(I) 

250 ISTSISH-NDIM^I 

kRITt{SCRC) {SUMUn, I-liNST0RE)*(SUM2n)', isl'iNPARl), 
, < 8 BIAS(I}»I»liMPARn 
NEWARC?,TRyB-;. 

RETURN 

ENTRX 6STIMB<SNREVjKN3TR, IKNSTR) 

DO 300 I«1,HdX2 

IUK = NUHR&R(SNREVU WK^STRi iKNSTR) 

IF<IUK.eQ,0} GO TO 300 
1,15KMSTR(UK*IKNSTRJ 
IF(U,GTrNPAt3I) GO TO 300 
DO 275 j?lij:npari 
I1=NPARI«I*U 
ii=iNDXM0nj:) + u 

275 SUNUI1)*SUM1U1)+1 lODO/SI ASSG(UJ«*»2 

300 CONTINUE 

310 ISTsiNDXNO(SSTART) 

I1=IST+NSTA8T 

CALL SXMINV^SUMl(Il), ')PC0 MiNPCOH,.D£LiTA(NSTART)) 

I2=1ST 

DO 350 IfNSTARTiNDiM 
DELTAU)f 0.9D0 
11=12+1 

DO 340 J^NSTART* I 

DELTA ( n ?DELTAn >+6UMl(U)«SUM2(U) 

340 ii=ii+NniM-y 

IFn,EQ,NDlPI) GO TO 348 
IPlsl+1 

DO 345 J»IP1;NDIM 
U-IST + J 

345 DELTAU)*DELTa(I )+SUMl(Il><»SUM2(J) 

348 BBlAS(n?8BUS(n+DELTA(n 
350 IST=IST+NDIM«I 
ENDFILE SCR0 
REWIND SCRC 
RETURN 

entry UPDAT0 

READtSCRC) tSUMKi ) »I»3#NSTORE)» (SUM?n),l*liNPARU» 
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^ 23 .. 

i 30 

; 3 i 

132 

133 


134 

139 

136 

137 
<MH> ^ H 


. (BBIAS(n;i 51 *NPARU 

11 = 0 ... . . 

00 5Q0 I?1,SPARI 
DEUTA( P*0.§D0 
. DO 490 LtNSTART^NDIM 
Ll = ni'«!U 

450 DELTAS U»DEt;’TA(U-SyMl(Ll)»I)EUTAXL>' 

BBlASU)fDECTA( n+B 81 AS(I). 

500 Il=U + NDl’M-l 
return 

END 

7 MEMORY EXRANDao, USE SLlMlTS OR CORE* OPTION FOR NEXT RUN 
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1 

2 

3 

4 

3 

6 

7 

8 

? 

10 

U 

la 

13 

11 

19 
16 

17 

18 
1? 

20 
21 
22 
23 


26 

27 

?8 

P.f 

30 

31 

32 

33 

34 

39 
38 

37 

38 
3V 

40 

41 

42 


CSYMINV 

C NAMR SYMINV 

C 

C PURPOSE -.TO. RECURblVElY FIND INKERSE OF SYMysTRlC MATRIX 

C . . 

C CALLING SeQ'JENcS C*ALL SY^lNVOUMUNyiMiNLlM’/DELTAJ 

c „ 

C SYMBOL T»PE description 

c 

c SUM! nP INPUT p uoi^&R rectangular Part of matrix 
C (1) TO R6 INVfcRTfcD 

C OUTPUT lower rectangular PART QP INVERTED MATRIX 

C 

C NDIM I input » dimension qf matrix 

c 

C NLIM I INPUT S dimension of partition TO BE INVERTED 

r* 

c delta op Scratch 

c <i) 

c 

c subroutines used UONE 

r* 

u 

C common SlOcnS none 

n 

C INPUT FILES NONE 

C 

C OUTPUT FILES NONE 

C 

C references 'GEODYN systems DEScRI^TIONi 

C VOLUME 1 « GEODYN DOCUMENTATION 

c 

c 


subroutine SyMtNV(SUMl,NDlM,NUM»DEL.TA) 
DOUBLE PRECISION SUMK 1 ) , DELTA< 1 > 

. C INITIALIZE RY FInDInG INVERSE OF 1 X 1 
SUMlCl)al,gD 0 / 3 UHl(l) 

IP(NUIM,EOtI) RETURN 


N 1 =NDIM -1 

C RFCURSIVELY F'lNQ INVERSE OF NXN KNOWI.NG 'INVERSE OF (N-UXlN-i) UNTIL 
C THE INVERSE OF AN NLIM X N^lM SQUARE PARTITION IS FOUND 
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43 


DO .400 N?2,HLIM 

44 


NHlsN-1 

45 


L1*0 

46 


DQ lOO L?liSiai 

47 


Jl = 0 

49 


DELTA(L)80 ,5D0 

4? 


DO 60 jsl,u 

50 


JL=J1+L 

51 


JN=J1+N 

52 


delta (U) =DEU:Ta (L USJMl ( Jl, )#SUM1 ( gN ) ■ 

53 

60 

gi=Jl + NDl^l~rf 

54 


IF(L.EQ,NMU GO To lOD 

55 


LPlsL*l 

56 


DO 00 JaLPi;NMl 

57 


JNeJl+N 

58 


JL«U*J 

59 


DELTA (U "DELTA a ) + SUMl(JU)<*SUMl(gN) 

60 

80 

vll»v»l + NOlM-U 

61 

100 

LlsLl+NDlM-U 

62 


J1=N 

63 


nn?ni*n 

64 


DO 150 g?l,SMl 

65 


SUMllNN)?SUHl(NN)-D£LTA<g)*SUMl{'gi) 

66 

150 

JlaJl+NOlM-y 

67 


SUM1(NN)»1.9D0/SUm1(NV) 

6§ 


Jl"N 

6? 


DO 200 J*1,MM1 

70 


SUMI(J1)?-DElTA( J)«SUM1(NN) 

71 

200 

Jl = Jl+N0:M-kj 

70 


I1 = N 

73 


DO 300 I?1,Sm1 

74 


JlsI 

75 


DO 250 

76 


SUHKJD^SUMK J1)«SUMICI1)»DELTAU) 

77 

250 

Jl"Jt+NDlM-y 

78 

300 

Il"U+NniM-l; 

7? 

400 

Nl=Nt+NOlM-S 

80 


RETURN 

81 


END 
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